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Abstract. We assess the extent to which the flux-limited IRAS PSCz redshift survey encapsulates the complete 
or major share of matter inhomogeneities responsible for the external tidal forces affecting the peculiar velocity 
flow within the Local Supercluster and its immediate surroundings. We here investigate this issue on the basis of 
artificially constructed galaxy catalogs. Two large unconstrained iV-body simulations of cosmic structure formation 
in two different cosmological scenarios form the basis of this study. From these iV-body simulations a set of galaxy 
mock catalogs is selected. From these a variety of datasets is selected imitating the observational conditions of 
either the local volume-limited Local Supercluster mimicking NBG catalog or the deeper magnitude-limited PSCz 
catalog. The mildly nonlinear dynamics in the "mock" Local Supercluster and PSCz velocities are analyzed by 
means of the Least Action Principle technique in its highly optimized implementation of the Fast Action Method. 
By comparing the velocities in these reconstructions with the "true" velocities of the corresponding galaxy mock 
catalogs we assess the extent and nature of the external tidal influence on the Local Supercluster volume. We 
find that the dynamics in the inner 30/i _1 Mpc volume is strongly affected by the external forces. Most of the 
external forces can be traced back to a depth of no more than 100/i _1 Mpc. This is concluded from the fact that 
the FAM reconstructions of the 100/i _1 Mpc PSCz volume appear to have included most gravitational influences. 
In addition, we demonstrate that for all considered cosmological models the bulk flow and shear components of 
the tidal velocity field generated by the external distribution of PSCz galaxies provides sufficient information for 
representing the full external tidal force field. 
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1. Introduction 

Migration flows of cosmic matter are one of the major 
physical manifestations accompanying the emergence and 
growth of structure out of the virtually homogeneous pri- 
mordial Universe. The cosmic flows displace matter to- 
wards the regions where ever more matter accumulates, 
ultimately condensing into the objects and structures we 
observe in the Universe. 

Within the gravitational instability scenario of struc- 
ture formation, the displacements are the result of the cu- 
mulative gravitational force exerted by the inhomogeneous 
spatial matter distribution of continuously growing den- 
sity surpluses and deficits throughout the Universe. This 
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establishes a direct causal link between gravitational force 
and the corresponding peculiar velocities. Given a suitably 
accurate measurement of peculiar velocities within a well- 
defined "internal" region of space, Vint-, we may invert 
these velocities and relate them to the inducing gravita- 
tional force. Hence, the source of the measured motions 
may be traced and possibly even reconstructed. In princi- 
ple, it may even allow us to infer the total amount of mass 
involved and thus estimate the value of the cosmological 
density parameter U, m and other fundamental cosmologi- 
cal parameters. 

The practical execution of such studies of cosmic ve- 
locity flows is ridden by various complicating factors. One 
major complication is that the cosmic regions in which pe- 
culiar velocities have been determined to sufficient accu- 
racy may have a substantially smaller size than what may 
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be deemed appropriate for a dynamically representative 
volume. Ideally, in order to account for almost the com- 
plete flow in our local cosmic neighbourhood we should 
have probed the density field in a sufficiently large cosmic 
volume. This should involve a region of space substantially 
superseding that of the characteristic scale of the largest 
coherent structures in the Universe. Only then the mag- 
nitude of the gravitational influence of inhomogeneities 
at larger distances will represent a negligible contribution 
and, as well, start to even out against each other. 

The size of this dynamically effective volume depends 
sensitively on the structure formation scenario which is 
prevailing in our Universe. Hence, it will be closely af- 
filiated to the spatial distribution, characteristic size and 
coherence scale of cosmic structures, and its size will there- 
fore be in the order of the scale of the largest pronounced 
structures in the Universe. Within the conventional struc- 
ture formation models, based on Gaussian initial den- 
sity and velocity fields, this is fully specified through the 
scenario's fluctuations power spectrum P(k). When the 
power spectrum involves a substantial large-scale compo- 
nent and the survey volume is rather limited we have to 
be aware of significant external influences. Although not 
yet exactly determined, observational evidence suggests its 
size to be in the range of « 100 — 200/i _1 Mpc (where h de- 
notes the Hubble constant in units of 100 km s _1 Mpc _1 ). 

An equally important consideration concerns the spa- 
tial resolution at which the velocity field is studied. 
Available samples of galaxy peculiar velocities extend out 
to reasonable depth of « 60ft. _1 Mpc. Yet, they involve a 
rather coarsely and inaccurately sampled cosmic velocity 
field. By absence of precise distance estimators, more accu- 
rately and densely sampled velocity information is there- 
fore mostly confined to a rather limited region in and 
around the Local Supercluster [LS]. As a consequence, 
most analyses of large-scale cosmic flows are necessarily 
confined to spatial scales at which the evolving cosmic 
structures are still residing in a linear phase of develop- 
ment. The dynamics in more advanced stages of cosmic 
structure formation are as yet poorly constrained by mea- 
surements. 

1.1. Cosmic Force Fields and Supercluster Dynamics 

In this work we wish to extend the analysis of cosmic 
flows to the more advanced evolutionary stages pertain- 
ing within supercluster regions. Only within the local cos- 
mic neighbourhood of our Local Supercluster, the quality, 
quantity and spatial coverage of the peculiar velocity data 
are sufficiently good to warrant an assessment of the cos- 
mic velocity field and the corresponding dynamics at a 
sufficiently high spatial resolution. On these quasi-linear 
or mildly nonlinear scales we hope to find traces of the on- 
set towards the more advanced stages of cosmic structure 
formation. In order for this to yield a meaningful and suc- 
cessful analysis, two major questions have to be addressed. 
Both form the main focus of this contribution. 



The first issue, that of the rather restricted sample 
volume, constitutes the major incentive behind this work. 
The volume of the galaxy catalog that best samples our 
local cosmic neighbourhood, the Nearby Galaxy Catalog 
(Tully 1988, hereafter [NBG]), is certainly substantially 
smaller than what may be considered dynamically repre- 
sentative. Any analysis of the (internal) velocity field in 
our Local cosmic neighbourhood should therefore take into 
account the impact of external gravitational influences. 

We focus on two related problems. In the first place, 
there is the need to quantify the effect of neglecting the 
external gravitational influence g ext when modeling the 
cosmic velocity field on scales comparable to that of the 
Local Supercluster. Various studies have attempted to de- 
termine cosmological parameters on the basis of a com- 
parison between modelled versus measured velocity field 
in the Local Universe (Tonry & Davis 1981, Tully & Shaya 
1984, Shaya, Peebles & Tully 1995, Tonry et al. 2000). 
For this it is crucial to understand in how far local den- 
sity perturbations may account for the local peculiar grav- 
ity field within the Local Supercluster. Directly related to 
this is the need to have a sufficiently accurate descrip- 
tion of the external force field g ex t, in terms of its nature 
and spatial extent, in order to properly model the total 
measured gravity field gtot- For studies intent on a com- 
parison of modelled versus observed peculiar velocities on 
scales larger than the Local Supercluster this is an essen- 
tial requirement (Faber & Burstein 1988, Han & Mould 
1990, Yahil et aZ.1991, Webster et al. 1997, Branchini et 
al. 1999). Similar considerations are equally relevant for 
the inverse problem, in which one attempts to infer the ex- 
ternal gravitational influence g ext from peculiar velocities 
measured within the Local Supercluster (e.g. Lilje, Yahil 
& Jones 1986, Lynden-Bell & Lahav 1988, Kaiser 1991, 
Hoffman et al. 2001). Indeed, both problems concerning 
the restricted galaxy sample volume have figured promi- 
nently in previous cosmic velocity field studies and were 
addressed in a variety of publications. However, usually 
these tend to discard the fact that the local cosmic region 
in which we have access to high quality velocity data has 
already reached an advanced quasi-linear dynamical state. 

Referring to the latter, the second major issue con- 
cerns the innovative way in which we evaluate the dy- 
namical state of superclusters. These structures reside in 
a mildly nonlinear evolutionary stadium, having evolved 
significantly beyond their initial linear phase. Unlike the 
vast majority of previous studies, we seek to probe into 
the more detailed and informative kinematic aspects of 
these structures. A conventional linear analysis will not 
be able to provide an adequate description, and for the 
most evolved circumstances not even the Zel'dovich ap- 
proximation (Zel'dovich 1970) may be expected to do so. 
In order to be able to optimally exploit the available ve- 
locity information - without suffering the loss of valu- 
able high-resolution information through a filtering pro- 
cess - we apply the Least Action Principle [LAP] formal- 
ism (Peebles 1989) for dealing with the individual galaxy 
velocities. To that end, an optimal implementation devel- 
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oped by Nusser & Branchini (2000, hereafter [NB]), the 
Fast Action Minimization [FAM] proved an essential tool. 

Elaborating on the first issue of the external gravita- 
tional influence, one of the as yet undecided issues is the 
extent to which a LAP analysis of a cosmological self- 
gravitating system is dependent on a proper representa- 
tion of the external gravitational influence. Various strate- 
gies have been followed, ranging from a complete neglect 
of external forces (Peebles 1989), or taking account of the 
influence of merely a few nearby objects (Peebles 1989, 
1990, Dunn & Laflamme 1993, 2001), towards methods 
involving approximate descriptions of external influences. 
The latter mostly incorporated the wider external influ- 
ence through a frozen, linearly evolving, external tidal field 
estimated on the basis of the present-day locations of an 
extended sample of objects deemed representative for the 
external matter distribution (e.g. Shaya, Peebles & Tully 
1995, Schmoldt & Saha 1998, Sharpe et al. 2001). This 
study does include the influence of force fields, but does 
so in a fully systematic and self-consistent fashion, enabled 
by the FAM method to take into account the evolution of 
the full sample of external matter concentrations. 

The principal conclusion of our study is that the grav- 
itational forces exacted by the matter inhomogeneities 
encapsulated by the IRAS-PSCz redshift survey sample 
(Saunders et al. 2000) are indeed able to account for all 
motions within our local Universe. In addition, we demon- 
strate that its external influence may almost exclusively 
be ascribed to the bulk and shear flow components. 

1.2. Strategy 

This study is based on a number of artificial galaxy sam- 
ples mimicking the properties of genuine catalogs. They 
consist of several well-defined and well-selected model cat- 
alogs of galaxies and galaxy peculiar velocities. These 
mock samples are extracted from a set of extensive N- 
body simulations: for the nearby Universe models they ad- 
here to the characteristics of the Nearby Galaxy Catalog, 
for the deep galaxy redshift samples they are modelled af- 
ter the IRAS PSCz catalog. These model samples allow us 
to thoroughly investigate the various strategies forwarded 
for a successful and conclusive analysis. 

Two sets of realistic mock catalogs of galaxies are ex- 
tracted. The first set, the "local" one, is meant to mimic 
the mass distribution within the LS as traced by galaxies 
in the Nearby Galaxy Catalog of Tully (1988). It consists 
of a volume-limited galaxy sample within a (spherical) 
interior region with radius 30/i _1 Mpc. Each of these in- 
terior samples is embedded within a larger mock sample, 
the "extended" sample which covers a larger cosmic re- 
gion. In addition to the interior volume-limited sample in 
the inner 30/i _1 Mpc they contain an enclosing outer flux- 
limited sample covering the surrounding spherical region 
located between 30/i _1 Mpc < x < 100/i _1 Mpc. This ex- 
terior sample mimics a flux-limited galaxy catalog whose 
characteristics are modelled after the IRAS PSCz sample. 



For both the "local" and "extended" mock samples we 
model the peculiar velocity field at the positions of the 
particles in the "local" cosmic region, i.e. for the objects 
out to a radius x < 30/i _1 Mpc. These model predictions 
result from the application of Fast Action Minimization 
method (Nusser & Branchini 2000). As our FAM recon- 
struction procedure only takes into account the gravi- 
tational forces between the particles in the mock sam- 
ples - i.e. does not include contributions from outside 
objects - the differences in predicted velocities between 
the "local" and "extended" samples will reflect the influ- 
ence of the mass concentrations in the surrounding region 
30/i _1 Mpc < r < 100/i _1 Mpc. The comparison with the 
corresponding TV-body velocities, representing the "real" 
velocities, will inform us in how far "sky-covering" samples 
of galaxies with a depth of « 100/i _1 Mpc may be expected 
to represent a proper cosmic region as far as its dynam- 
ics are concerned. The strategy of analyzing and compar- 
ing the velocity models obtained from the small "local" 
mock catalogs, the large "extended" mock catalogs, and 
the "real" velocities in the original A-body simulations, 
will yield a solid understanding of the effect of neglecting 
the externally induced peculiar gravitational acceleration 
g ext (see Eqn. (1)). 

When analyzing a dataset of galaxy peculiar veloci- 
ties in the local Universe. The analysis of the larger PSCz 
mimicking catalogs should elucidate if and to what ex- 
tent such samples will be able to account for g ext . If the 
galaxies in these samples indeed appear to be responsible 
for the major share of the external forces, we may feel re- 
assured to use the PSCz sample of galaxies for a proper 
representation of g e xt- 

One aspect of this question concerns the investigation 
of the question whether the external tidal influence may 
be explicitly framed in an analytical approximation con- 
sisting of a dipolar and quadrupolar term. Our fully self- 
consistent FAM reconstructions, in which the "extended" 
mock catalogs are processed with the inclusion of all ex- 
ternal matter concentrations, enable us to estimate the 
bulk and shear components in the induced "local" galaxy 
motions. By comparing the resulting velocity fields in the 
"local" and "extended" samples we will be able to judge 
the quality of the approximate methods, and quantify 
and investigate the possible presence of systematic trends 
throughout the "local" cosmos. 

To account for possible systematic effects due to global 
cosmology, the mock galaxy catalogs are extracted from 
A-body simulations in two different cosmic structure for- 
mation scenarios. One involves a KCDM Universe with 
a characteristic large-scale dominated power spectrum, 
while the other concerns a tCDM cosmology. The more 
small-scale dominated character of the latter leads to a 
different character of its gravitational field fluctuations, 
the smaller coherence scale of the density field fluctua- 
tions yielding a comparatively smaller influence of the ex- 
ternal (quadrupolar) tidal field (the induced bulk flows 
are similar, as the smaller tCDM fluctuations are exactly 
compensated by the larger mass involved). The resulting 
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comparisons of FAM velocity field reconstructions are ex- 
pected to reflect these velocity field differences. 

In the end, this study of artificial galaxy samples 
should allow us to appreciate the manifestations of the real 
physical effects we wish to grasp. In this, we also should 
learn how to deal with the complications due to the host of 
measurement uncertainties which beset the observational 
data. The scope is to quantify the systematic errors which 
might have affected similar, local, comparisons based on 
real data and to judge whether the information on the 
external mass distribution available to these analyses is 
indeed sufficient to account for g ex t- 

1.3. Outline 

In the next section, we will elaborate on the astrophysi- 
cal background of this study, the study of velocity flows 
on cosmological scales, and in particular the issue of in- 
ternal and external gravitational influences. Ensuingly, we 
address the specific problem of treating the dynamics and 
related cosmic motions within mildly nonlinear structures 
such as the Local Supercluster. This brings us to a brief 
exposition on the LAP analysis for dealing with the com- 
plications of mildly nonlinear orbits and the technical is- 
sue of the FAM technique which allows us to apply this 
to a system composed of many objects. Special emphasis 
is put on the inclusion of external gravitational influences 
within the LAP /FAM formalism. In section 4 we describe 
the cosmological setting of the simulations on which this 
study has been based. As a guidance towards interpret- 
ing our calculations, we address a variety of theoretical 
aspects and predictions concerning cosmic velocity fields 
in these cosmological scenarios. The basis of this work 
is the set of two "parent" N-body simulations and the 
mock catalogs extracted from these simulations, forty in 
total. They are presented in section 5. In the subsequent 
sections we present the results obtained from the various 
FAM computations. In section 6 we analyze the veloc- 
ity vector maps for the FAM reconstructions. These maps 
allow a direct and visually illuminating appreciation of 
the effects wc: wish to address. This is followed by a first 
quantitative assessment in section 7. This consists of a 
comparison between the FAM velocity field reconstruc- 
tions of the Local Supercluster volumc(r < 30ft. -1 Mpc), 
the FAM reconstructions for the corresponding PSCz sam- 
ple and the complete "real world" N-body velocity field. 
The comparison is mainly based on a point-by-point evalu- 
ation through scatter diagrams of velocity-related quanti- 
ties. To encapsulate these results into a spatially coherent 
description of the large scale external velocity and grav- 
ity field, in section 8 we turn to a decomposition of the 
peculiar velocity field into multipolar components. In par- 
ticular, we demonstrate that a restriction to its dipolar 
and quadrupolar components, i.e. the bulk flow and veloc- 
ity shear, does represent a good description. Thus having 
looked at the issue of cosmic velocity fields from different 
angles, the summary of section 9 will focus on the reper- 



cussions of our analysis and its relation to the study of 
the (relatively nearby) surrounding matter distribution. 
On the basis of these conclusions we provide a descrip- 
tion of the various projects which follow up on this work, 
together with some suggestions for additional future work. 

2. Cosmic Flows: 

probes of cosmic matter distribution 

2.1. the Large-Scale Universe: linear flows 

Over the past two decades a major effort has been di- 
rected towards compiling large samples of galaxy pecu- 
liar velocities. These samples made it possible to obtain 
a rather impressive view of cosmic dynamics on scales 
> 10/i _1 Mpc. In particular the Mark III catalog, with 
an effective depth « 60/i _1 Mpc, stands as a landmark 
achievement (Willick et al. 1997a, also see Dekel 1994 and 
Strauss & Willick 1995). Further progress has been en- 
abled by the assembly of additional and partially comple- 
mentary samples of galaxy peculiar velocities, of which the 
SFI late-type galaxy and ENEAR early-type galaxy sam- 
ples are noteworthy examples. The SFI Catalog of Peculiar 
Velocities of Galaxies (Giovanclli et al. 1997a, Giovanelli 
et al. 1997b, Haynes et al. 1999a and Haynes et al. 1999b) 
consists of around 1300 spiral galaxies with I-band Tully- 
Fisher (TF) distances, out to cz < 7500 km s -1 . The 
ENEAR sample (da Costa et al. 2000) is an equivalent 
sample of around 1600 early-type galaxies, out to a dis- 
tance cz < 7000 km s _1 , with D n — a distance estimates 
available for nearly all of them. Tracing cosmic motions 
over larger volumes of space is a rather more cumbersome 
affair and attempts to do so are mainly based upon the 
peculiar motions of galaxy clusters. The claim of a puz- 
zlingly large flow over scales of 150/i _1 Mpc by Lauer & 
Postman 1994 could not be corroborated. Nonetheless, 
flows on such large scales may indeed be a reality, as 
has been inferred from the far better defined "Streaming 
Motions of Abell Clusters" (SMAC) sample of Hudson 
et al. 2001. They did recover a bulk flow in the order of 
687 ± 203 km s _1 , of which 225 km s _1 may arise from 
sources at a distance larger than 100/i~ 1 Mpc (Hudson 
et al. 2004). One prime objective of most analyses of 
these large samples of peculiar velocities has been the de- 
termination of the cosmological mass density parameter 
n m (Davis, Nusser & Willick 1996, Willick et al. 1997b, 
Willick & Strauss 1998, Nusser et al. 2000, Branchini et 
al. 2001). Such assessments are based on a comparison of 
observed velocities to a model velocity field. A basic re- 
quirement for obtaining self-consistent estimates of fi m is 
that the velocity samples concern a "representative" vol- 
ume of space. 

However, even while such studies appear to succeed 
in attuning the large-scale matter distributions and veloc- 
ity fields in a reasonably self-consistent fashion, doubts 
remain with respect to a variety of practical and sys- 
tematic problems. Firstly, in these comparisons the ran- 
dom errors on the observed velocities are substantial, 
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much larger than those in the structure formation models. 
Considerable effort has been directed towards quantifying 
and minimizing errors on the observed peculiar velocities 
(e.g. Dekel 1994, Strauss & Willick 1995, and references 
therein). These involve random measurement errors as 
well as more subtle systematic, yet reasonably well under- 
stood, errors. Secondly, there remain various systematic 
effects which have not been addressed and corrected for 
in an equally convincing fashion. Even though they also 
tend to play a role with respect to the model predictions 
they are often overlooked. 

A major systematic factor concerns the incomplete in- 
formation on the spatial mass distribution within the re- 
gion of the sample itself. This prevents an adequate treat- 
ment of artifacts due to the incomplete sky coverage and 
limited depth of the available samples, and effects system- 
atic errors stemming from luminosity and density effects. 
These systematic errors are usually accounted for by using 
large, all-sky redshift surveys, such as the Optical Rcdshift 
Survey of Santiago et al. (1995) or the 1.2 Jy and PSCz 
surveys of IRAS galaxies (Fisher et al. 1995, Saunders et 
al. 2000). In particular when using IRAS based surveys 
the effects of incomplete sky coverage are greatly reduced. 

Even more problematic for a successful handling of 
luminosity and density related effects is our incomplete 
knowledge with respect to the relationship between the 
observable galaxy distribution and the underlying mass 
distribution. By absence of a compelling theory of galaxy 
formation this "galaxy bias" is usually encapsulated in 
heuristic formulations. The rather ad-hoc and possibly un- 
realistic or inadequate nature of the latter may seriously 
affect the significance of the inferred conclusions. Most 
studies make the simplifying assumption of a galaxy pop- 
ulation fairly tracing the underlying density field. This is 
usually embodied in a global and linear "galaxy bias" fac- 
tor. A large variety of results suggest that this may be 
a reasonable approximation on scales in excess of a few 
Megaparsec. Moreover, while this bias may be problem- 
atic in the case of early-type galaxies, it has proved to 
be quite successful with respect to the later type galaxies 
which figure prominently in IRAS based samples (Verde 
et al. 2002). 

2.2. Internal and External influences 

Unlike most studies of cosmic flows which seek to assess 
and analyze the nature and source of dynamical influ- 
ences within a confined region of space, we try to get 
an impression of the cosmic dynamics on mildly nonlin- 
ear scales of only a few Megaparsec. We focus on the 
Local Supercluster region and its immediate neighbour- 
hood. The galaxy sample of the NBG catalog is taken to be 
representative for this region. Because the catalog entails 
a volume which is substantially smaller than what may be 
considered dynamically representative, the peculiar veloc- 
ities of the galaxies are partially due to the gravitational 
action by outside matter concentrations. That is, the pe- 



culiar velocities are not only due to the gravitational force 
induced by the matter concentrations within the "inter- 
nal" survey volume V, n t , but also reflect the gravitational 
influence by the "external" matter density distribution, 
g ext ■ Because it does not constitute a truly representative 
volume of the Universe, a meaningful dynamic analysis of 
the Local Universe on the basis of the NBG sample is sub- 
stantially complicated by its limited depth, which is one 
of the major systematic problems besetting the analysis of 
virtually all available surveys of galaxy peculiar velocities. 
Theoretical models of peculiar velocities nearly always in- 
volve the implicit assumption of the mass being homoge- 
neously distributed outside Vi nt , so that its gravitational 
effect may be neglected. Even in the case of having a suffi- 
ciently large volume at one's disposal, this approximation 
is only valid in the central part of Vi n t , certainly not near 
its edges. 

The distinction between external versus internal grav- 
itational force may be best appreciated by noting that 
the total (peculiar) gravity field gtot(x) is the netto sum 
of the individual contributions by all patches of matter 
throughout the visible Universe. At any position within 
the internal volume Vi nt , we may then decompose the full 
gravitational field into an "internally" induced component 
gint and an "externally" generated contribution g ex t, 

gtot(x) = gmt(x) + gext(x) . (1) 

In this way we have defined the internal gravitational force 
gint as the integrated contribution from the density fluc- 
tuations <5(x) within the volume Vint, while the external 
gravitational force g ext concerns the combined gravita- 
tional force generated by the density fluctuations outside 
the realm of Vi nt , so that 

. , mH 2 r , . r , . . (x' - x) 

g tot (x,t) = — — / dx'«5(x',i)^ 1 + 

3ttH 2 f , , c . , , (x' - x) 

/ dx'J x'.t > (2) 

8?F JVext l X - X l 

The peculiar velocities of galaxies within V in t bear the 
mark of both the acceleration due to the matter concen- 
trations within the volume itself, gint, as well as that of 
the combined gravitational influence of the external mass 
distribution, g ex t- A comparison of predicted internally 
induced velocities with the observed local velocity field 
should therefore enable us to infer the magnitude and na- 
ture of the external field gext- This analysis is usually 
facilitated by the fact that the fine details of the exter- 
nal force contribution are largely negligible. The contri- 
butions by the various external matter concentrations to 
the combined gravitational force mostly average out such 
that what remains noticeable is mainly confined to the 
low order components of the multipole decomposition of 
gext- This can be most readily appreciated from a descrip- 
tion of the external gravitational force field in terms of its 
successive multipole components. When we expand g ex t 
around some central location in Vint ~ here defined to be 
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the origin of the coordinates x - we find that to second 
order 

3 

9ext,i(x) = 9bulk,i - ^ %j Xj . (3) 
j=l 

In this, we assume that the additional divergence term 
j(V • gext)xi has been embedded into the (zcroth) order 
monopole term. In essence, it corresponds to a "breathing 
mode" affecting the "local" Hubble expansion within the 
volume, and therefore can not possibly be inferred from 
the local measurement of the internal gravity field gi„ t . 

The leading term in the overall external gravitational 
acceleration is the bulk gravity term gtmlk,i- This dipole 
term constitutes the uniform acceleration of the matter 
within V int , 

f , x' 

Sbulk = — ^dx<Kx,i )RF (4) 

Evidently, when considering peculiar velocities relative to 
the centre of mass inside the volume Vint instead of ab- 
solute velocities this constant vector disappears. The first 
term whose magnitude and configuration is independent 
of the reference frame is the quadrupolar term %j . 

If the contribution to the (peculiar) gravitational po- 
tential by the external mass inhomogeneities is <j> ex t, the 
quadrupolar tidal tensor %j is the trace-free part of 
d 2 (j) ex t/dxidxj, evaluated at the centre of Vi. It is deter- 
mined by the external matter distribution through 

% 3 (t) = — dx 5 (x , t) | — b 1 (5) 

The integral expressions for the dipole and quadrupolc 
components of the external gravity field (eqn. 4 and 5), 
illustrate that it is unfeasible to exploit the observed lo- 
cal cosmic velocity field to recover the detailed and com- 
plete spatial distribution of the external matter inhomo- 
geneities. On the other hand, it does indicate how it is 
that we can infer some overall characteristics of the ex- 
ternal matter distribution from an analysis of the local 
velocity field. From this we may extract interesting and 
significant information on the nature and even distribu- 
tion of the large scale cosmic matter distribution and set 
constraints on the values of some of the fundamental cos- 
mological parameters. The pioneering work by Lilje, Yahil 
& Jones (1986) in which the velocity field of the Local 
Supercluster was exploited to infer the presence of a ma- 
jor external source of gravitational attraction has shown 
the potential of this approach. Ultimately, it inspired the 
analysis of Lynden-Bell et al. (1988) that lead to the dis- 
covery of the Great Attractor. 

3. Cosmic Flows: 

the mildly nonlinear dynamics of superclusters 

Even though a structure's evolution may have progressed 
to a dynamical stage at which a first-order description of 



cosmic velocity fields will no longer be adequate, it may 
still be possible to find a direct link to the structure's ini- 
tial configuration. This is in particularly true for the early 
and mildly nonlinear phases of evolution. The exemplary 
archetype of a structure in which such mildly nonlinear 
circumstances are prevalent is that of superclusters, the 
filamentary or wall-shaped elements of the cosmic foam- 
like matter distribution. 

Over the past two decades intriguing foamlike patterns 
have gained prominence as a prime characteristic of the 
cosmic matter distribution. The first indications for the 
actual existence of a foamlike galaxy distribution were 
provided by CfA2 redshift slices (de Lapparent, Geller & 
Huchra 1986) and established as a universal cosmic pat- 
tern with the Las Campanas redshift survey (Shectman et 
al. 1996). With the arrival of the large recent and ongoing 
systematic galaxy redshift surveys, the 2dF galaxy red- 
shift survey (« 250, 000 redshifts, Colless et al. 2003, also 
see e.g. Colless 2004 and Tegmark et al. 2002 for a discus- 
sion on clustering in the 2dFGRS) and the Sloan Digital 
Sky Survey (SDSS, will determine « 1,000,000 redshifts, 
see e.g. Zehavi et al. 2002 and Tegmark et. al. 2004 for 
an overview of present-day status wrt. galaxy clustering), 
we may hope to have entered the stage in which we will 
be enabled to explore the formation and the dynamics of 
these characteristic spatial structures in the cosmic mat- 
ter distribution. The typical elements of the cosmic foam - 
filamentary and wall-shaped superclusters - are precisely 
at the youthful yet mildly nonlinear phase of development 
mentioned earlier. They were identified as such within the 
context of Zel'dovich' "pancake" theory of cosmic struc- 
ture formation (see e.g. Shandarin & Zel'dovich 1989). The 
significance of the cosmic foamlike network for the un- 
derstanding of the process of cosmic structure formation 
has since been generally recognized. This may be appreci- 
ated from the widespread use of the concept of the 'cosmic 
web', coined by Bond, Kofman & Pogosyan (1996) in their 
study of the dynamics underlying its formation (see Van 
de Weygaert 2002 for a recent general review). 

Mildly nonlinear cosmic features such as superclusters 
have recently turned their initial co-expansion into a gen- 
uine physical contraction (or are on the brink of doing so) , 
marking the emerging structure as a truely identifiable en- 
tity. Once it has surpassed this "turn-around" stadium the 
complexity of its internal kinematics quickly increases. At 
first moderately, and ultimately dramatically as the viri- 
alization process advances, the matter orbits inside the 
structure become more and more complex. Even in the 
more moderate early phases of this process, an appropri- 
ately sophisticated treatment of the mildly nonlinear dy- 
namics appears to be a necessary requisite for any study 
based upon kinematic information. In and around emerg- 
ing nonlinear structures a simple linear analysis for recon- 
structing initial conditions will therefore no longer suffice. 
In other words, a sufficiently detailed and profound un- 
derstanding of the emergence of these key elements in the 
cosmic matter distribution cannot be obtained without 
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the development of a more elaborate technique for the 
analysis of cosmic velocity fields. 

3.1. Structure formation: mildly nonlinear dynamics 

A linear analysis simplifies the dynamical evolution of a 
system into an initial conditions problem. It implies the 
reconstruction of the primordial density and velocity field 
by means of a simple linear inversion of the observed mat- 
ter distribution and galaxy peculiar velocity field. Such an 
approach may even be followed towards a slightly more 
advanced stage. The Zel'dovich formalism, a Lagrangian 
first-order approximation for gravitationally evolving sys- 
tems, has been remarkably successful in describing the 
early nonlinear evolution of a supercluster (for a review, 
see Shandarin & Zel'dovich 1989). Substantially surpass- 
ing its formal linear limitations, it proved to be a highly 
versatile medium for analyzing and explaining the overall 
spatial morphology and characteristics of emerging struc- 
tures. The Zel'dovich approximation elucidated and ex- 
plained qualitatively the fundamental tendencies of gravi- 
tational contraction in an evolving cosmos. Perhaps most 
noteworthy this concerned the tendency of gravitational 
collapse to proceed anisotropically, together with its pre- 
dictive power with respect to location and timescales of 
the first phase of collapse into planar mass concentrations, 
"pancakes". This offered the basic explanation for the 
foamlike morphology of the cosmic matter distribution, 
stressing its existence many years in advance of its discov- 
ery through observational programs to map the galaxy 
distribution (for an extensive review of various nonlin- 
ear approximation schemes seeking to expand upon the 
Zel'dovich approximation see Sahni & Coles 1995). 

In line with the above, the Zel'dovich approximation 
proved a highly versatile tool for the analysis of the cos- 
mic matter flows. It was successfully applied to the nonlin- 
ear situation of mixed boundary conditions - tested and 
calibrated using iV-body simulations - by Nusser et al. 
(1991) and Nusser & Dckel (1992). However, its validity 
remains restricted to the early stages of nonlinearity at 
which there is still a linear and direct relation between ve- 
locity and gravity field. Once matter inside the emerging 
structures starts to reach densities so high that local inter- 
actions become dominant, the Zel'dovich scheme quickly 
ceases to lose its applicability. Once matter elements start 
to cross each each others path the interaction between the 
nonlinear matter concentrations within the realm of the 
contracting structure will more and more deflect the or- 
bits away from their initial linear trajectory. The linear 
kinematics of the Zel'dovich approximation will therefore 
no longer be able to follow the orbits of the matter ele- 
ments. Higher order approximations based on perturba- 
tion theory have been advocated in order to follow such 
more advanced nonlinear circumstances. However, the im- 
provement over simple first order Zel'dovich approxima- 
tion turns out to be quite limited and not warranting the 
effort invested at each successive perturbation step. This 



is particularly so as with the onset of nonlinearity the rate 
at which successive perturbative orders terms become sig- 
nificant rapidly accelerates. 

3.2. Least Action Principle in Cosmology 

In more advanced nonlinear circumstances we encounter 
a more generic dynamical situation than a simple initial 
value problem. Typically, one seeks to compute the ve- 
locity field consistent with an observed density structure 
at the present epoch or, reversely, one deduces the den- 
sity from the measured peculiar galaxy velocities. In the 
case of generic systems, the dynamical evolution repre- 
sents a mixed boundary condition problem. This implies 
the system to be sufficiently constrained by complement- 
ing the incomplete dynamical information regarding the 
initial conditions with that pertaining to the dynamical 
state of the system at a different epoch. While iV-body 
codes are particularly concerned with the ideal circum- 
stances usually embodied in terms of initial value prob- 
lems, a different kind of technique needs to be invoked to 
exploit the typical mixed boundary information yielded by 
observations. 

A more profound and direct exploitation of the avail- 
able information to follow the physics of such a cosmo- 
logical nonlinear system was forwarded by Peebles (1989, 
1990). He pointed out that finding the orbits that satisfy 
initial homogeneity - and by implication vanishing initial 
peculiar velocities - and match the (present-day) observed 
distribution of mass tracers constitutes a mixed-boundary 
value problem. Such problems lend themselves naturally 
to an application of Hamilton's principle. This naturally 
leads to the formulation of the Least Action Principle (also 
known as "Numerical Action Method"), based on a varia- 
tional analysis of the action S of an isolated system of M 
particles, which at a cosmic expansion factor a(t) is given 

by 

S = jf Ldt = jT dt \\^ 2 ^ - mi (j>{Ti) } , (6) 

in which L is the Lagrangian for the orbits of particles with 
masses nii and comoving coordinates Xi and correspond- 
ing peculiar gravitational potential 0(x). For a system of 
particles interacting by gravity alone, embedded within 
a uniform cosmological background of density pb(t), this 
yields the following explicit expression for the action S, 

s = L * I?— br) + 7gR^r 

2 -I 

+ -irGpta^m^ 2 ^ (7) 

The exact equations of motion for the particles are then 
obtained from finding the stationary trajectories amongst 
the variations of the action S subject to fixed boundary 
conditions at both the initial and final time. 

Confining oneself to a feasible approximate evaluation 
in this Least Action Principle approach, one describes the 
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Fig. 1. 2D projected re- 
constructed FAM orbits 
for different values of Nf . 
The black dots represent 
the final (present) po- 
sitions for each object. 
The solid lines indicate 
the trajectories followed 
by the objects as a func- 
tion of time. The top- 
left panel shows FAM re- 
constructed orbits with 
Nf = 1 (Zel'dovich ap- 
proximation) . The top- 
right with Nf = 2, the 
lower-left Nf = 3 and the 
lower-right Nf = 6. 



orbits of particles, Xj(t), as a linear combination of suit- 
ably chosen universal functions of time with unknown co- 
efficients specific to each particle presently located at a po- 
sition x^o- For instance, by using the linear growth mode 
D(t) as time variable (Giavalisco et al. 1993, Nusser & 
Branchini 2000), one can parametrize the orbit x^(D) of 
a particle as 

71=1 

where the functions q n {D) form a set of Nf time- 
dependent basis functions. The factors Cj n are then a 
set of free parameters, whose value is determined from 
evaluating the stationary variations of the action. 

The functions q n {D) satisfy both two orbital con- 
straints, necessary and sufficient to define solutions 
in agreement with evolution in the context of the 
Gravitational Instability theory for the formation of struc- 
ture in the Universe: q n (l) = ensures that at the present 
time the galaxies are located at their observed positions 
Xj(l) = Xj :0 and hiriD-^ D s ' 2 q n (D)0(D ) = guaran- 



tees vanishing peculiar velocities at early epochs which, 
in turns, ensures initial homogeneity. 

3.3. Fast Action Minimization 

The successful application of the Least Action Principle 
towards probing the kinematics and dynamics of an evolv- 
ing cosmological system depends to a large extent on the 
specific implementation. This will be dictated by the char- 
acteristics of the physical system. In order to enable a 
meaningful LAP analysis of large samples of galaxies, like 
the Local Universe samples studied in this work, an op- 
timized procedure is necessary for dealing with the large 
number of objects. Nusser & Branchini (2000) developed 
an optimized version of Peebles' original LAP formalism, 
the Fast Action Minimization method. The various opti- 
mization aspects of the FAM implementation proved to be 
crucial for our purposes. The relevant optimization hinges 
on four major aspects of the FAM scheme. 

The first FAM improvement involves the choice of time 
basis functions q n (D). Its convenient choice of time ba- 
sis functions yields a simple expression for the action of 
the system and for its derivatives with respect to Cj )Tl . 
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Both quantities relate to the internal gravity term gj„ t of 
the system. Once the action and its derivatives are evalu- 
ated numerically the minimum of the action is determined 
by means of the conjugate gradient method (Press et al. 
1992). The orbits of the system are then fully specified 
through the set of parameters found in correspon- 
dence to the minimum. 

Closely related to the first aspect is that of tuning 
the choice of the time basis functions q n {D) such that 
only a limited number Nf of basis functions is needed to 
successfully parameterize the orbits of the system. This is 
in particular beneficial to the the physical configuration we 
are studying here, involving Megaparsec scale dynamics 
characterized by quasi-linear or mildly nonlinear motions. 

Note that using the growth factor D as time variable 
makes the equations of motions almost independent of the 
value of il m (Nusser and Colberg 1998). As a consequence 
FAM orbits and peculiar velocities in a generic £l m uni- 
verse can be obtained by appropriate scaling those assum- 
ing a flat cosmology. 

A final major aspect of the FAM implementation 
involves the efficient computation of the internal (self- 
consistent) gravity gi nt from the particle distribution in 
the sample. To this end, the gravitational forces act- 
ing on the particles at the different epochs are com- 
puted by means of the TREECODE technique (Bouchet & 
Hernquist 1988). By proceeding in this fashion, the FAM 
method is able to reconstruct the orbits of « 10 4 — 10 5 
mass tracing objects back in time. This makes FAM nu- 
merically fast enough to perform a large number of orbit 
reconstructions, essential for performing the intended sta- 
tistical analysis presented in the following sections. 

In this work we use Nf — 6 basis functions to pa- 
rameterize the orbits, choosing a tolerance parameter 
tol = 1CP 4 to search for the minimum of the action S and 
setting a softening parameter of 0.27/i _1 Mpc to smooth 
the gravitational force in the TREECODE. Orbit search- 
ing in dynamically relaxed systems is a difficult exercise 
since one has to choose among the many solutions found 
at the extrema of the action. However, since the purpose 
of FAM is to investigate large scale dynamics dominated 
by coherent flows rather than virial motions, our evalua- 
tions translates into an orbit search restricted to solutions 
which do not deviate too much from the Hubble flow i.e. to 
the simplest orbits that represents the minima of the ac- 
tion. Therefore, we set the initial guess for Ci <n according 
to linear theory prescription and search for the minimum 
of the action to avoid multiple solutions found a station- 
ary points which typically describe passing orbits (Peebles 
1994). We have checked that this choice of parameters is 
optimal in the sense that decreasing tol, increasing Nf or 
changing the input set of Cj jTJ does not modify the final 
results appreciably. 

Distortions in the resulting FAM-predicted peculiar 
velocities mainly arise from two systematic artifacts 
(Branchini, Eldar & Nusser 2002). One is the discrete sam- 
pling of the mass distribution within Vi nt ■ The second, and 
overriding one, is the failure of FAM in reproducing the 



virial motions within high-density regions that is a direct 
consequence of having considered solutions that represents 
perturbations to the Hubble flow. This deficiency of the 
FAM reconstructions is clearly illustrated by the resid- 
ual velocity vector maps (see eq. 20) in Fig. 7 and Fig. 8 
(bottom row). These show the velocity vector differences 
between the "real" measured, i.e. A-body, velocities and 
the corresponding FAM reconstructions (here based on 
cither the galaxy distribution in a 30ft- _1 Mpc central re- 
gion or the extended 100/i _1 Mpc region). The maps show 
how the largest residuals are the ones found in the high 
density regions: although the FAM30 and FAM100 veloc- 
ity fields do show pronounced velocities near these regions 
they are not the proper "real" virialized velocities they 
should have been. The residual fields thus underline the 
fact that FAM's inaptitude to deal with regions character- 
ized by large virial motions. Instead, in those situations it 
will lead to a false prediction of coherent inward stream- 
ing velocities, an effect pointed out by Nusser & Branchini 
(2000) and which can be also noted in our images when 
carefully studying them. 

Finally, for practical reasons, since we are merely in- 
terested in measuring the effect of external gravity fields 
we make a further simplifying hypotheses. We ignore red- 
shift distortion effects by working in real space In this 
respect, we should point out that extensions of the action 
principle method allowing a direct processing of redshift 
space information have been proposed and shown to work 
(Phelps 2000, Phelps 2002, also see Sharpe et al. 2001). 

3.4. The role of biasing 

In this work we perform orbit reconstructions by assum- 
ing that all the mass of the systems is associated to point 
mass objects. More explicitly, we are making two different 
hypotheses. The first one is that we are able to identify a 
set of objects that trace the underlying mass density field 
in an unbiased way. The second one is that the internal 
structure of these objects is irrelevant for our reconstruc- 
tion purposes. 

The first assumption hardly applies to real galaxies 
that are most likely to be biased tracers of the mass dis- 
tribution, as indicated by the relative bias between galax- 
ies with different luminosities, colors and morphological 
type (Loveday et al. 1995). However, if galaxies and mass 
particles share the same velocity field so that the biasing 
relation remains constant along the streamlines, then the 
problem can be easily circumvented by specifying the bi- 
asing scheme at the present epoch (Nusser and Branchini 
2000). 

Within the standard lore of galaxies embedded in a 
virialized halo of dark matter that grow through hierar- 
chical merging of smaller systems, neglecting the internal 
structure of objects is an assumption that is best justified 
a posteriori by showing how well Numerical Action meth- 
ods can reproduce N-body velocities. Although the good- 
ness of this assumption has been quantified by a number 



10 



E. Romano-Diaz et al.: External Tidal Impact on Local Supercluster Dynamics 



of numerical tests (e.g. Nusser and Branchini (2000) and 
Branchini Eldar & Nusser (2002)) little effort has been 
devoted to understand why numerical action methods can 
accurately reconstruct the velocity field of a large N-body 
simulation. 

One of the reason for this success is that only ~ 5% of 
the points used in our reconstructions, that were randomly 
selected from the N-body simulation, belong to virialized 
regions where FAM reconstruction fails. Fortunately, the 
locality of this "virial effect" allows us to partially circum- 
vent this problem by applying a modest spherical tophat 
smoothing of 2/i _1 Mpc to the FAM-predicted velocities. 
This tophat filter has been specifically important for the 
quantitative aspects of our study, where such systematic 
problems may sort distorting conclusions. This smooth- 
ing has been invoked in quantitative comparisons between 
FAM and iV-body velocities presented in this work, in par- 
ticular in the regression analyses. 

Little is known about the ability of numerical action 
methods to reconstruct the orbits of virialized systems. 
Indeed, when applied to extended objects rather than 
point masses, numerical action methods follow a single 
center of mass point per virialized objects, completely 
neglecting its merging history. Some argument can be 
given to back our choice of neglecting the internal struc- 
ture of virialized objects. First of all, after tracing back 
the merging history of virialized halos in N-body exper- 
iments a simple visual inspection reveals that particles 
ending up in the same halo at z=0 are contained within 
regions with simple boundaries at high redshifts. As a 
consequence, high order terms in the gravity potential 
about the halo center of mass are probably rather small. 
This probably minimize the role of major mergers whose 
rate for galaxy-size halos peaks in the redshift range 2-4 
(Volonteri, Haardt and Madau 2003) while peculiar mo- 
tions mostly develop at z < 2 (Branchini and Carlberg 
1994). These qualitative arguments clearly need to be con- 
firmed by appropriate numerical analyses similar to that 
of Branchini and Carlberg (1994) but extending out to 
scales of cosmological interest. 

3.5. Ordered reconstructions 

To obtain an idea of the level of improvement obtained 
through the use of successively higher order FAM evalua- 
tions, Figure 1 depicts 2D projections of the corresponding 
FAM particle orbit reconstructions within a local spheri- 
cal volume of « 30/i _1 Mpc. The black dots indicate the 
positions for each object in the sample, while the lines em- 
anating from each dot represent the computed trajectories 
followed by these objects as they moved towards their 
present location. The illustrated configuration is taken 
from one of constructed mock catalogs, and resembles that 
of the Local Universe (see section 5.2.1). Each successive 
FAM reconstruction is based on the same (present-day) 
particle distribution. The four frames correspond to suc- 
cessively higher order FAM approximations, involving an 



increasing number Nf of basis functions q n (D). The top- 
left panel shows FAM reconstructed orbits with Nf = 1, 
which in fact corresponds to the conventional first or- 
der Zel'dovich approximation and thus represent the or- 
bits that would have been obtained by the PIZA method 
(Croft & Gaztahaga 1997). These are followed by panels 
with Nf = 2 (top right), Nf = 3 (bottom left) and Nf = 6 
(bottom right). They show a clear and steady improve- 
ment towards the Nf — 6 FAM evaluation. Testing proved 
that even higher order computations did not yield im- 
provements significant enough to warrant the extra com- 
putational effort. 

In summary, the galaxy orbits in our Local Universe 
environment are found at a minimum of the action which 
is not too far, yet different, from linear theory predictions. 
The FAM technique thus yields a significant modification 
of the recovered galaxy orbits and peculiar velocities for 
configurations that evolved well beyond the linear regime 
(see e.g. Figure 1). Potentially its ability to deal with non- 
linear circumstances might even prove of benefit to recover 
sets of cosmological initial conditions satisfying nonlinear 
observational constraints at the present day, which indeed 
has been suggested by Goldberg & Spcrgel 2000. 

3.6. LAP and External forces 

The original cosmological Least Action Principle formula- 
tion by Peebles ( 1989) considered a fully self-consistent, 
i.e. isolated, system of point masses. For practical rea- 
sons, the original implementation had to be restricted to 
systems of at most a few dozen objects. Almost exclu- 
sively, the Local Group of galaxies formed the focus of 
these LAP studies (Peebles 1989, 1990, 1994, Dunn & 
Laflamme 1993). 

While these studies did indeed yield a substantial 
amount of new insight into the dynamical evolution of 
the Local Group, the issue of incorporating the dynami- 
cal influence exerted by external mass concentrations re- 
mained a major unsettled question. External forces do 
represent a significant component of the dynamics of the 
Local Group, as had been shown by Raychaudhury & 
Lyndcn-Bcll ( 1989). They established beyond doubt that 
the Local Group cannot be considered a tidally isolated 
entity, and demonstrated that the Local Group is acted 
upon by an appreciable quadrupolar tidal force. The re- 
sulting tidal torque appears to be responsible for the large 
angular momentum of the Local Group as a whole, as 
Dunn & Laflamme (1993) showed in an elegant and pio- 
neering analysis using orbits computed by the LAP vari- 
ational method. They confirmed that the tidal influence 
of the external matter distribution is indeed essential to 
explain its present angular momentum. 

In the course of time various strategies emerged to in- 
clude external dynamical influences. The nature of these 
methods are mainly set by the character of the physical 
system under consideration, and to some extent was de- 
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pendent on the available computational resources. Three 
strategies are outlined below. 



3.6.1. Directly Including External Masses 

To incorporate the external tidal influence within the 
LAP analysis the work by Peebles (1989, 1990, 1994), 
Peebles et al. 2001 and Dunn & Laflamme (1993) at- 
tempted to identify a few principal external mass concen- 
trations which would be responsible for the major share 
of the external gravitational influence. While in his first 
LAP study Peebles (1989) considered the Local Group 
mainly as an isolated system, sequel studies (Peebles 1990, 
Peebles et al. 2001) attempted to assess the possible ex- 
ternal influence by neighbouring matter concentrations. 
In Peebles 1990 he attempted to condense the exter- 
nal tidal force into two nearby mass concentrations, the 
Sculptor and Maffei group, each modeled as a single mass. 
Both were incorporated as 2 extra particles, with prop- 
erly scaled masses, within the action S in order to take 
them along in a fully self-consistent variational evaluation. 
A similar approach was followed by Dunn & Laflamme 
(1993), be it that they included five galaxies/groups in the 
local cosmic neighbourhood which arguably contribute a 
significant torque on the Local Group. Also in a later ap- 
plication (Peebles et al. 2001) this approach was followed, 
be it with an extensive outer region between 4/i 75 1 Mpc and 
20/i 75 1 Mpc whose mass distribution was condensed into a 
coarse sample of some 14 major external objects. 

This "self-consistent" strategy is feasible to pursue 
within the context of the original, computationally inten- 
sive, LAP implementation. This approach may therefore 
be followed in LG resembling situations in which a few 
objects suffice to represent the main aspects of a system's 
dynamical evolution. On the other hand, cosmic systems 
of a considerably larger scale than the Local Group would 
in general be too demanding for. Supercluster sized re- 
gions, with scales of up to a few tens of Megaparsec, count 
many more individual objects than a galaxy group. These 
systems have also not yet reached a formation stage so 
advanced that they have already largely decoupled from 
the global Hubble expansion, so the resulting external 
gravitational influence is usually the shared responsibil- 
ity of a large number of external matter concentrations. 
Accounting for the large-scale tidal field would quickly 
become prohibitively expensive in terms of the computa- 
tional effort for conventional LAP analyses. 



3.6.2. Inserting External Tidal Potential 

An alternative strategy is to incorporate the external 
gravity in the LAP scheme via an approximate expres- 
sion for the external contribution. This may be most di- 
rectly achieved by inserting an extra external tidal poten- 
tial term $udai(t) in the action S (eqn. 6). As on suffi- 
ciently large, Megaparsec, scales we may expect this term 



to evolve according to linear gravitational instability per- 
turbation growth, 



®tidal(t) = D^t) <S>tidal(t ) 



D(t) 
a(t) 



$ tidal (to), 



(9) 



in which <&tidai{to) is the tidal term at the present cosmic 
epoch and = (D/a) the linear growth term for the 
gravitational potential (the growth factor and cosmic 
expansion factor a(t) are set to be equal to unity at the 
present epoch, D(to) = a(to) = 1.). Thus, instead of evolv- 
ing it self-consistently along with the considered system, 
the external field is determined at one epoch - usually the 
present one - and then incorporated as an extra linearly 
growing gravity field term &udai(t) in the action S: 

o f'" r \ ^ m i« 2 ( dxA 2 G m-^ mirrij 

s = I dt IE — (irj + 

+ -irGpba 2 ^™^ -$tidal(t)\ ■ (10) 

There are various possibilities to compute the tidal poten- 
tial term $udai, usually from the current mass distribu- 
tion. One option is to compute it directly from a sample 
of M ext external objects which is deemed responsible and 
representative for the major share of the external tidal 
force field, 
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Note that none of these external objects (j = 1, • • • , M ex t) 
is taken into account as far as the action of the system and 
the computation of their orbits is concerned, except for 
their "passive" role in determining $ 'tidal- An alternative 
approach is to insert an approximate analytical expres- 
sion for Qtidai , in particular one including the dipolar and 
quadrupolar contributions, d and T, to the tidal potential, 



(to) = E d 
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+ - Xj 
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T X,; 



(12) 



Equivalently, one may chose to insert the corresponding 
expressions directly into the expression for the derivative 
of the action with respect to an expansion coefficient, 
dS/dC° ! n = 0, evidently equal to zero within this vari- 
ational approach. 

The first, "direct", procedure (eqn. 11) was followed 
by Shaya, Peebles & Tully (1995), who for the pur- 
pose of studying the velocity field within the surrounding 
30/i _1 Mpc modelled the relevant external mass distribu- 
tion after the distribution of rich Abell clusters from Lauer 
& Postman (1994). To some extent, Sharpe et al. (2001) 
operated along the same lines, be it that they added the 
resulting tidal term directly to the reconstructed veloci- 
ties produced by the LAP procedure. However, while in 
principle exact, such a concentrated and static mass dis- 
tribution may involve considerable uncertainties and can 
be highly sensitive to the uncertainties in the location of 
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a few dominant point masses. As this spatial point distri- 
bution is supposed to form a suitable model for the un- 
derlying large scale matter distribution this may be even 
more worrisome. 

Potentially more elegant may therefore be the mod- 
elling of a smooth tidal field along the line of the second 
procedure (eqn. 12), as suggested by Schmoldt & Saha 
(1998). The corresponding dipolar and quadrupolar term 
may then be based on the best available determinations of 
these parameters. On the other hand, when the LAP vol- 
ume is comparatively large, the analytical approximation 
may represent an oversimplification of the force field, ne- 
glecting potentially important local variations within the 
external force field. 

3.6.3. Selfconsistent and Direct FAM approach 

The indirect "potential" approach which we described 
above (eqn. 11 or eqn. 12) may not properly account for 
the temporal evolution of the external field in the case 
of nonlinear ly evolving systems. The formalism assumes 
a static, merely linearly evolving, gravitational potential. 
However, the matter concentrations which generate the 
external tidal forces will themselves get displaced as the 
cosmos evolves. These displacements may be relatively mi- 
nor for distant masses, but for the more nearby entities 
this may be entirely different. A detailed treatment of the 
external mass distribution will be necessary when the in- 
fluence of the nearby external objects on the evolution of 
small "interior" regions is comparable to or even dominant 
over the selfgravity of the region. It will be equally cru- 
cial to follow the detailed whereabouts of nearby matter 
concentrations in the case of a large "interior" region in 
which a marked contrast between the central regions and 
the outer realms may result in a significantly different dy- 
namic evolution. 

This prompted us to follow an alternative and direct 
approach, a fully self-consistent strategy in which also 
the external matter concentrations are accounted for in 
the computation of the system of evolving particle orbits. 
Alongside that in the "local" region for which we seek to 
reconstruct the velocity field, also the system of objects 
in the exterior regions (30/i _1 Mpc < r < 100/i _1 Mpc) are 
considered. Non- uniform manifestations of the external in- 
fluence can only be included by pursuing such a direct and 
systematic approach. It is only through the availability of 
the FAM technology that we were enabled to do so for 
a Megaparsec system consisting of a large number of ob- 
jects. 

4. Cosmological scenarios 

The mock catalogs on which we apply our Fast Action 
Minimization analysis are extracted from iV-body simula- 
tions in two different cosmological settings. Their charac- 
teristics, in terms of their relevant parameters, are listed in 
Table 1. The table also lists the simulation specifications. 
The first scenario is a flat ACDM model with a cosmolog- 



ical constant term f2A,o = 0.7 (Clo = 0.3, Oo + ^A,o = 
1.0, T = 0.25, n = 1). The second model is a tCDM 
Einstein-de Sitter (fi = 1.0, A ,o = 0,T = 0.25, n = 1) 
model, motivated by the decaying particle model proposed 
by Bond & Efstathiou (1991). Both scenarios were cho- 
sen to be viable with respect to the current observational 
constraints, implying similarities in many overall proper- 
ties and appearances, yet with some significant differences 
with respect to their dynamical repercussions. This may 
provide indications on whether the galaxy motions in our 
local cosmic neigbourhood do contain information on the 
structure formation scenario. 

In both cases the amplitude of density fluctuations is 
normalized on the basis of the observed abundance of rich 
galaxy clusters in the local universe. This abundance de- 
pends on the magnitude of the matter field fluctuations 
on the mass scale characteristic for galaxy clusters. This 
translates into a dependence on the amplitude of den- 
sity fluctuations on cluster scales modulated by the mean 
global matter density. A variety of studies (e.g. White, 
Efstathiou & Frenk 1993, also see Eke, Cole & Frenk 1996) 
found that in order to yield the present-day cluster abun- 
dance the amplitude of density fluctuations in spheres of 
radius 8/i _1 Mpc, as, and fio are related by 

cr 8 = 0.551V ' 6 , (13) 

The resulting power spectra are depicted in Figure 2 (top 
left). On all scales, the density fluctuations in the rCDM 
scenario, represented by the dotted lines (for both P(k) 
(green lines) and k 3 P(k) (blue lines)), are less pronounced 
than those of the ACDM scenario: the two power spec- 
tra have a similar shape and differ by a simple scaling 
factor over the entire wavelength range. Visually, this is 
immediately reflected in the stark differences between the 
spatial galaxy distribution in the resulting mock catalogs. 
Figure 3 provides such a visual comparison. It shows the 
"external" PSOz catalog mimicking galaxy distribution in 
three mutually perpendicular central slices in the case of 
the ACDM scenario (top row), together with the same set 
of frames for a rCDM mock galaxy catalog (bottom row). 
On all scales, the rCDM galaxy distribution looks con- 
siderably more uniform than that in the ACDM Universe. 
Not only is the clustering of galaxies in the ACDM sce- 
nario more pronounced, it also delineates considerably 
larger structures, a manifestation of the power spectrum's 
amplitude at the corresponding large wavelengths. 

Because the higher average matter density in the f^o = 
1 tCDM Universe does almost fully compensate for the 
lower amplitude of the density fluctuations the resulting 
gravity and velocity perturbation fields in the ACDM and 
tCDM scenarios are very similar. The velocity power spec- 
tra k 3 P v (k) are shown in the bottom lefthand panel of 
Figure 2: their functional dependence is the same over 
the entire wavelength range. The larger mass correspond- 
ing to a given density excess in the rCDM Universe evi- 
dently effects a stronger gravitational force. The resulting 
large scale motions scale as f(flo) K ^o' 6 - This happens 
to be almost exactly the inverse of the average density 
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Table 1. iV-body simulation parameters. Column 1: cosmological model. Column 2: flo mass density parameter. Column 3: A, 
cosmological constant parameter. Column 4: F, power spectrum shape parameter. Column 5: <r&, density perturbation amplitude 
spectral normalization. Column 6: size of the computational box in /i _1 Mpc. Column 7: number of particles in the simulations. 

Cosmology fio A T as boxlength N bj 

ACDM 0.3 0.7 0.25 1.13 345.6 192 3 

tCDM 1.0 0.0 0.25 0.55 345.6 192 3 



perturbation amplitude scaling (cq. 13), which is propor- 
tional to fijT 0,6 (eqn. 13). While this is exactly the factor 
involved in the normalization of the power spectrum, in 
terms of as, the lower level of density fluctuations gets pre- 
cisely compensated by the higher amount of mass involved 
with them. This can be directly observed from the velocity 
power spectra P v (k) for the two scenarios (Fig. 2, lower 
lefthand frame). The velocity power spectra for both sce- 
narios are exactly equal over the entire wavelength range, 
both in functional dependence as well as in amplitude. 
Note that also the gravity perturbations in the ACDM sce- 
nario are substantially stronger than those in the rCDM 
cosmology: because they scale with \£IHq and the am- 
plitude of the density perturbations, which according to 
eq. 13 is oc O ' 6 , the average peculiar gravitational accel- 
eration is proportional to Vt <}A H^. 

The comparison between k 3 P(k) (Fig. 2, top panel, 
blue lines) and k 3 P v (k) (Fig. 2, bottom panel) in the same 
figure shows the shift of the velocity perturbations, with 
respect to the density perturbations, towards a more large- 
scale dominated behaviour. This follows directly from the 
continuity equation, connecting the velocity and density 
perturbations such that the velocity power spectrum re- 
lates to P(k) through P v (k) oc P(k) k~ 2 . 

The large-scale behaviour of the (linear) velocity per- 
turbation field immediately illuminates the difficulty in 
tracing the full array of matter inhomogeneities responsi- 
ble for the cosmic motions within a specific cosmic region. 
To account for all noticeable contributions it is necessary 
to probe out to large depth. This is manifestly evident for 
the first order component in the externally induced flow, 
the "bulk flow" v bu i k . A measure for the expected bulk 
flow within a (tophat) spherical region of size Rth, 

v bulk (x.) ee / dx'v(x')W T H(x-x',i? TH ), (14) 
Jv 

is represented by the (root square) average value a v , whose 
value may be inferred from the Fourier integral 



= H f(n ) 



(15) 




dk 



3 P(k)W^(k)fc-2 



In these relations, Wth (x, Rth) and WthQl) are the ex- 
pressions for the tophat window filter, spatially and in 



Fourier space, and er_i is the spectral moment aj for 
j = — 1 (see Bardeen et al. 1986, henceforth BBKS). 

How substantial the large scale origin of the bulk flow 
is may be readily appreciated from figure 2 (top centre). 
Because the linear character of fluctuations on large scales, 
the spectral a v (eq. 16) does provide a reasonable order- 
of-magnitude estimate of the magnitude of the large-scale 
bulk motions. The figure shows the estimated bulk flow 
amplitudes, « (J y . CIS <X function of the (tophat) window 
radius: the bulk flow is clearly a large scale phenomenon, 
converging only very slowly towards large spatial scales. In 
both the ACDM scenario and the tCDM scenario the ex- 
ternally induced bulk flow on a scale of 30/i _1 Mpc will be 
in the order of 200 km s _1 . Of this overall bulk flow, more 
than 100 km s _1 has to be ascribed to inhomogeneities on 
scales exceeding 200/i _1 Mpc ! When assessing the motions 
in a local volume of 30/i _1 Mpc radius, in terms of relative 
external contributions, inhomogeneities on a scale larger 
than 100/i _1 Mpc still contribute more than 25% of the 
total while the ones larger than 200/i _1 Mpc are still re- 
sponsible for more than 10% (see Fig. 2, centre bottom). 
We should therefore expect to find substantial external 
contributions in the ACDM and rCDM simulations. Note 
that this relative contribution to the bulk flow, the "bulk 
convergence", is defined as the relative contribution by 
matter perturbations within a radius Rth to the exter- 
nally induced bulk flow on a scale of 30/i~ 1 Mpc (the size 
of the NBC volume): 



ulk 



1 - 



<?v(Rth) 
(30/i _1 Mpc) 



(16) 



The second order aspect of the velocity field which we 
seek to study is the induced velocity shear Sij , 



1 { dvi dvj 1 1 ._ . . 

2W + ai -3 (v ' v) ^' 



(17) 



Also the velocity shear reveals interesting and distinguish- 
ing differences between the ACDM and the rCDM sce- 
nario. In the linear regime the expected magnitude of the 
shear tensor Sij, on a tophat scale Rth, may be evalu- 
ated through its direct proportionality to the tidal shear 




Fig. 2. Spectral characteristics for the two studied cosmological scenarios. Upper left: Power spectrum P(k) (green lines) and 
k 3 P(k) (blue lines), for the ACDM (continuous lines) and rCDM (dotted lines) cosmological scenarios used in this study. Lower 
left: velocity power k 3 P v (k), the amount of velocity power per logarithmic waveband interval. Upper centre: expected bulk flow 
velocities ( km s _ ) over a spherical (tophat) volume of radius Rth ■ The solid line is the real theoretical prediction, the dashed 
one concerns the restricted power spectrum represented in the simulation box. Centre low: the convergence of the bulk flow, 
defined according to eqn. 16. Upper right: Expected shear flow ( km s _1 Mpc _1 ) over a spherical (tophat) volume of radius Rth- 
Similar to bulk flow, the solid line concerns the full theoretical power spectrum while the dashed corresponds to the restricted 
power spectrum represented in the simulation box. Lower right: convergence of the shear flow, defined according to eqn. 16. 
Note that for the righthand frames the corresponding ordinate axis are the ones on the righthand side of the plotbox, while for 
the central frames it is the one on its lefthand side. 



Tij . Quantifying Sij by means of its (root square ) average 
a s (van de Weygaert & Bertschinger 1996), we find 



a s = H f{fl )a {R TH )^^- (18) 

in which the (dimcnsionless) spectral parameter 7 = 
af/ao 02 is defined through the th , 1 st and 2 nd spectral 
moments Uj (see BBKS 1986). The predictions for the 
two cosmological scenarios are shown in topright frame of 
Fig. 2. With respect to the bulk flow there is a marked 
difference in coherence scale: the major contributors to 
the tidal shear are located at considerably closer distances 
than the sources of the bulk flow (fig. 2, cf. lower right 
with lower centre). Most of the shear inducing matter in- 
homogeneities are found within a radius of 100/i -1 Mpc, 
accounting for more than 95% of its value (fig. 2, lower 
right). On a scale of 30ft. - 1 Mpc we expect an external tidal 
shear of « 7 km s^Mpc -1 for both the ACDM scenario 
and the rCDM model. 



5. Mock catalog Construction and Analysis 

5.1. N -body simulations 

The two A-body simulations used in this work were car- 
ried out by Cole et al. (1998) within the context of an ex- 
tensive study of PSCz catalogue resembling galaxy mock 
samples in a large variety of cosmological structure forma- 
tion scenarios. They consist of 192 3 particles in a computa- 
tional box of 345.6 ft. Mpc. They are dynamically evolved 
using an AP3M code in which the force is smoothed with 
a softening parameter of 0.27 /i _1 Mpc. 

The purpose of this study is a demanding task for tru- 
ely representative A^-body simulations. The A^-body sim- 
ulations should provide an optimal compromise between 
a high mass resolution on the small scale side and, on 
the large-scale side, a cosmic volume large enough to be 
dynamically representative. The large dynamic range re- 
quirement involves a mass resolution refined enough to re- 
solve mass entities comparable to galaxies. This translates 
into an average inter-particle separation that needs to be 
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Fig. 3. Central slices through a mock galaxy catalog realization. Depicted is the galaxy distribution in the external (r > 
30 h~ 1 Mpc) volume. For each catalog, we show the x — y (left), x — z (centre) and y — z (right) plane. Top row: ACDM mock 
catalog; Bottom row: rCDM catalog. 



smaller or comparable to that of galaxies in real observa- 
tional catalogues. On the other hand, the simulations have 
to extend over a cosmic volume which is large enough to 
incorporate the major share of the gravitational influence 
exerted by the inhomogeneous cosmic matter distribution. 
Given the slow convergence of the bulk flow and its large 
coherence scale this is particularly challenging, and will be 
in the order of several hundreds of Megaparsec (see discus- 
sion in the previous section and fig. 2). Although hardly 
any current TV-body simulations would fully fulfill the dy- 
namic range requirements, the used TV-body simulations 
do appear sufficiently adequate for a meaningful investi- 
gation of the relevant systematic trends and effects. This 
remains true in a qualitative sense, even though on the 
basis theoretical arguments (see e.g. Fig. 2) and observa- 
tional indications (e.g. Hudson et al. 2004) we know there 
may be substantial bulk flow contributions stemming from 
even larger spatial scales. 



In this respect it is important to note is that the mock 
catalog realizations in this work are constrained by the 
finite size of the simulation box. The practical repercus- 
sions of being confined to a limited simulation volume may 
be inferred from the dashed curves in Fig. 2. They show 
the corrections to the expected bulk flow and velocity 
shear predictions (solid curves) when only the inhomo- 
geneities in the restricted volume of the 345.6 ft, _1 Mpc 
simulation box are incorporated. Because perturbations 
on scales exceeding the fundamental scale of the box are 
absent, the realized power spectrum has a rather sharp 
and artificial large-scale cutoff: the limited boxsize Lt ox 
implies a cutoff in the power spectrum at low wavenum- 
ber kbox = 2ft/Lbox- From Fig. 2 we can conclude that 
this correction is particularly apt for the bulk flow, pre- 
dictions for the velocity shear seem hardly affected. As a 
consequence, on scales over w 100 h~ 1 Mpc the bulk flows 
in the realized TV-body simulations will be severely re- 
pressed and far from representative. Although large-scale 
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mode adding procedures have been proposed to partially 
remedy this situation (Tormen & Bertschinger 1996 and 
Cole 1997), our rCDM and ACDM simulations did not 
include such MAP (mode adding procedure) extensions. 
Conclusions with respect to the convergence of the FAM 
reconstructed velocity flows should therefore be referred 
to with respect to the suppressed velocity power spec- 
trum indigenous to our A^-body simulations (notice that 
this dynamic range issue is truely cumbersome to nearly 
any study attempting to assess velocity flows in computer 
simulations). 

5.2. Mock Catalog Construction 

From the full A^-body simulations we extract mock cata- 
logs made to resemble the local Universe. The ACDM and 
rCDM A^-body simulations are processed through spec- 
ified observational masks to imprint the required char- 
acteristics on the resulting mock catalogs. We distinguish 
two types of mock catalogs. From each A^-body simulation 
we extract ten different "local" mock catalogs mimicking 
the NBC catalogue and, with these "local" samples rep- 
resenting their interior, ten different "extended" samples 
resembling the PSCz catalogue. 

The "local" class of mock samples is meant to sample 
the mass distribution within a 30 h~ 1 Mpc region in and 
immediately around the Local Supercluster. These cata- 
logs constitute volume-limited galaxy samples mimicking 
the Nearby Galaxy Catalog of Tully (1988). Mock catalogs 
of the second type are designed to account for the mass 
distribution out to distances of 100 /i _1 Mpc. These "ex- 
tended" samples represent flux-limited samples, for which 
we take the IRAS PSCz galaxy catalog (Saunders et al. 
2000) as template. The PSCz sample is not only ideal for 
our purposes in that it covers one of the largest volumes 
of the Universe amongst the available galaxy redshift sur- 
veys, but also in that it concerns a survey covering a large 
fraction of the sky and involves a well-defined uniformity 
of selection. Assuming that on large linear scales IRAS 
galaxies define an unbiased tracer of the underlying dark 
matter, Hamilton ctal 2000 found that its real-space power 
spectrum is consistent with that of a COBE-normalizcd, 
untilted, flat ACDM model with Cl m = 0.3 and J7a = 0.7. 

In both flux-limited and volume limited samples the 
mass of mock galaxies have been rescaled to the value 
of TO = of the parent N-body simulations, listed in 
Table 1. 

In Table 2 we have listed the main characteristics of 
all mock catalogs used in this work. . 

In constructing the mock samples galaxies were identi- 
fied with N-body particles selected randomly, exclusively 
according to the catalog selection criteria. Therefore, we 
did not attempt to include bias descriptions to model pos- 
sibly relevant differences in the spatial distribution of dark 
matter and galaxies. This is different from the original use 
of the simulations (Cole et al. 1998), in which various bias 
descriptions were invoked to construct artificial galaxy 



Table 2. Characteristics of the mock catalogs. Column 1: label 
of the catalog. Column 2: cosmological model. Column 3: for 
each mock catalog type: number of generated catalogs. Column 
4: average number of particles (N ob j) in each catalog. Column 
5: external radius of the catalog in units of h~ 1 Mpc. 



Cosmology 


Set 


Number of 




Maximum 






catalogs 




radius 


ACDM 


NBC 


10 


2740 


30 




PSCz 


10 


10900 


100 


rCDM 


NBG 


10 


2803 


30 




PSCz 


10 


11207 


100 



samples whose two-point correlation function and large- 
scale power spectrum largely matched that of the APM 
survey (Maddox et al. 1996). The analysis of the small- 
scale nonlinear power spectrum of the PSCz by Hamilton 
& Tegmark ( 2002) even implies the bias on small scales 
to be very complex, involving a scale-dependent galaxy- 
to-mass bias. We, however, prefer not to include an extra 
level of modelling prescriptions. Our interest concerns the 
kinematics and dynamics of the matter distribution in the 
Local Universe, and the velocities of galaxies are thought 
to reflect these almost perfectly: they are mere probes 
moving along with the underlying dark matter flows, ir- 
respective of their particular bias relation with respect to 
the dark matter distribution. The sole strict assumption 
is therefore that of having no velocity bias (Carlberg et 
al. 1990), which on the large-scale Megaparsec scales at 
hand should be a more than reasonable approximation. 

5.2.1. Mock NBG catalogs 

The mock NBG catalogs are obtained by extracting spher- 
ical volumes of 30/i _1 Mpc from the A^-body simulation 
particle distribution. The positions of the spheres in the 
parent simulations are not random but chosen to mimic 
as close as possible the characteristics of a Local Group 
look-alike region. Therefore, each mock catalog is centered 
on a particle moving at a speed of 625 ± 25 km s -1 , re- 
siding in a region in which the shear within 5/i _1 Mpc is 
smaller than 200 km s _1 , where the fractional overdensity 
measured within the same region ranges between -0.2 and 
1.0. 

The velocity vector of the central particles defines 
a Galactic coordinate system and a Zone of Avoidance. 
Particles within the Zone of Avoidance are removed and 
substituted with a population of synthetic objects dis- 
tributed using a random-cloning technique (Branchini et 
al. 1999). The Zone of Avoidance [ZA] in the mock samples 
is designed to mimic that of the PSCz catalog (Saunders 
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et al. 2000) and is smaller than the one of the real NBG 
catalog (Tully 1988). 

Each spherical region contains on average 2 x 10 4 par- 
ticles. This set of particles is randomly resampled in order 
to produce an unbiased catalog of around 2800 objects, a 
number that matches that of the galaxies in the real NBG 
catalog (Table 2). This procedure preserves, within shot 
noise errors, density fluctuations and thus does not alter 
the orbit reconstruction. 

These NBG mimicking mock catalogs define volume- 
limited galaxy samples, so that the number of objects 
within a distance x therefore increases as x 3 . This is indeed 
what the resulting realizations yield, as may be discerned 
from the central part of the corresponding histogram in 
Figure 4 (x < 30ft _1 Mpc). 

5.2.2. Mock PSCz catalogs 

The second set of mock catalogs was obtained by carv- 
ing out spherical regions of radius 100ft _1 Mpc from the 
A^-body simulations. Each of these new mock samples is 
centered on the same central position as that of corre- 
sponding NBG mock catalogs, with which they share the 
objects within the central 30ft _1 Mpc. 

While the central 30ft" x Mpc region coincides with the 
NBG mock sample, the particle distribution in the exter- 
nal region (30ft _1 Mpc < x < 100ft _1 Mpc) is supposed to 
mimic that of galaxies in the flux-limited IRAS PSCz cata- 
log. To achieve this the objects beyond 30ft _1 Mpc were se- 
lected from the A^-body particle samples according to the 
PSCz selection function used by Branchini et al. (1999): 

/ x 2 \~ 13 

i>(x) = Ax- 2a M + —J ifx>x s . (19) 

In this expression x is the distance of the galaxy 
in ft _1 Mpc, while the parameters have the values a = 
0.53, (3 = 1.8, x s = 10.9/1" ^pc, and x* = S^^Mpc 
(Branchini 1999). The selection function ip(x) defines the 
relative number density of the galaxy sample with respect 
to the number density of the A^-body particle simulation. 
While the number density of PSCz galaxies is equal to 
that of the particles in the simulation at 10.9ft~ 1 Mpc, the 
amplitude A of ip(x) is normalized such that ^(10. 9) = 1. 

Two additional steps concern the treatment of "Zone of 
Avoidance" objects and the evening of the matter density 
throughout the full external sample volume. A first step is 
the processing of sampled objects in the Zone of Avoidance 
such that the resulting sample conforms to a reality resem- 
bling situation. The ZA "removal+substitution" is imple- 
mented in the same way as in the case of the NBG mock 
catalog construction, with the replacement achieved with 
the same random-cloning technique. Finally, in order to 
guarantee a uniform average mass density throughout the 
volume, the mass of the objects in the flux-limited external 
object sample (30ft _1 Mpc < x < lOOft^Mpc) has been 
scaled by the inverse of the selection function (j){x). 




Fig. 4. Object counts as a function of distance in the PSCz 
mock catalogs. Upper panel: the two histograms show the av- 
erage counts in the mock samples for the rCDM and ACDM 
cosmological models. The continuous line show the expected 
counts. Lower panel: the histograms show the fractional dif- 
ference between the observed and expected counts averaged 
over all PSCz mock samples in the two cosmological models 
explored. 

5.2.3. Mock catalog realizations 

From both figure 5 and figure 6 one can obtain an im- 
pression of the spatial context of the local NBG mock 
sample within the wider environment of the surrounding 
100ft _1 Mpc PSCz sample. To visually appreciate the se- 
lection criteria of the catalogs, and their interrelationship, 
figure 5 shows a three-dimensional view of one set of the 
ACDM mock catalogs, extracted from the particle distri- 
bution in the A^-body simulation of structure formation in 
a ACDM scenario. Emanating from the full PSCz + NBG 
mimicking galaxy samples in the top righthand cube is 
a row of two cubes showing the content of the external 
PSCz mock catalog (righthand) and the content of the 
central NBG mock catalog (lefthand) . Figure 6 elaborates 
on this, and shows the projected particle distribution in 
the same PSCz + NBG mock catalog (left panel) while 
focusing in on the central region (right panel). The circle 
(left panel) indicates the boundary of the volume-limited 
region comprised by the mock NBG galaxy sample, which 
in the right panel has been enlarged to show the corre- 
sponding velocity field within this NBG region. Velocities 
of objects within a 10ft _1 Mpc thick slice are shown by 
means of arrows whose size is proportional to the ampli- 
tude of the galaxy velocity components within this slice. 

As a matter of test, we checked the distance distribu- 
tion of the resulting mock galaxy samples. The histograms 
of the resulting mock catalog distributions are shown in 
Figure 4. The upper part of Figure 4 shows the num- 
ber of galaxies - averaged over all PSCz mock catalogs 
for both two cosmological models - as a function of dis- 




Fig. 5. Compilation of cubic subregions and various mock subsamples of the object distribution in one of the ACDM PSCz 
mock catalogs. The full mock catalog of objects, the PSCz look-alike volume and the NBG look-alike volume out to a radius of 
100/i _1 Mpc, is shown in the large lefthand cubic box (with size 200/i -1 Mpc). The boxes along the top row and those along the 
bottom row represent mutually exclusive object samples. The lower ones represents the inner cosmic region (x < 30/t _1 Mpc), 
the top row represent a corresponding part of the exterior PSCz region (30ft~ 1 Mpc < x < 100/i _1 Mpc). Emanating from the 
complete cubic volume and running along the top row are two successive boxes which zoom in on the external PSCz mock 
catalog objects, with cubic sizes of 160/i _1 Mpc (middle right cube) and 80/i -1 Mpc (lower right cube) respectively. The lower 
row of two cubic volumes the central (inner) NBG mimicking object sample. The full inner sample, comprising a sphere of 
radius 30/t -1 Mpc is shown in the lower central cubic volume, whose size of 80/i _1 Mpc is equal to the central PSCz box, while 
the righthand cube focuses in on the 40fo _1 Mpc size central region. 



tance x over the full range x < 100/i -1 Mpc. Clearly vis- 
ible is the discontinuity at 30/i _1 Mpc, marking the tran- 
sition from volume-limited NBG-like region to the PSCz- 
like flux-limited outer region. For comparison, the solid 
line shows the theoretically expected counts (Eq. 19). The 
generated mock samples appear to match the expected 
distance distribution rather well. This is further under- 



lined by the fractional difference between observed and 
expected counts, N b s /N cxp — 1, shown in the lower part 
of Figure 4. The fractional difference between mock sam- 
ples displays a perfect featureless scatter pattern: Poisson 
noise free of systematic effects. 
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Fig. 6. 2D projection of the particle distribution in one of the ACDM PSCz mock catalogs (left panel). The inner circle divides 
the NBG look-alike volume limited region from the PSCz flux limited one. The right panel show the 2D projected velocities in 
a slice of 10/i _1 Mpc cut through the inner part of the catalog. 



5.3. Mock Catalog Analysis 

5.3.1. FAM velocity field reconstructions 

On the basis of the FAM reconstructions of the galaxy 
velocities and the comparison with the true velocities - 
i.e. those in the original A^-body simulation - we assess to 
what extent the matter distribution within the confines 
of each different mock galaxy sample does contribute to 
the total velocity of the galaxies. In these idealized cir- 
cumstances of the ./V-body world the galaxy positions and 
velocities are known to perfect accuracy, thus circumvent- 
ing the need to investigate the effects of measurement er- 
rors and deceptive systematic biases in the galaxy peculiar 
velocities. This should provide us with a better under- 
standing of the nature and magnitude of genuine physical 
influences. 

Three different velocities are accorded to each galaxy 
located within the "local" spherically shaped 30/i~ 1 Mpc 
NBG region. The first velocity is that of the "true" , .TV- 
body velocity. For each of the in total 20 NBG mimicking 
galaxy mock catalogs, the FAM reconstructions produce 
two additional velocity estimates. One FAM velocity re- 
sults from the application of the FAM analysis to the re- 
stricted inner 30/i _1 Mpc NBG-like region itself. The sec- 
ond FAM based velocity is obtained on the basis of the 
FAM analysis on the extended, "full", 100/i _1 Mpc PSCz 
survey resembling sample (in which the "local" NBG sam- 
ple occupies the interior central region). In the following, 



we will indicate these FAM velocities by the names of 
FAM 30 and FAM100 velocities. 

5.3.2. FAM30 versus FAM 100 reconstructions 

The mutual comparison between each of the three differ- 
ent galaxy velocities - the FAM30, the FAM100 and the full 
./V-body velocities - is expected to yield abundant infor- 
mation on the dynamics and development of the structure 
in the interior 30/i -1 Mpc region: 

The FAM30 velocities are the galaxy velocities which 
would have been the product of the combined gravita- 
tional interaction of - solely - the matter concentrations 
within the central 30ft. -1 Mpc volume. Any deficiency with 
respect to the "real" ./V-body velocity of each galaxy has 
to be ascribed to the gravitational impact of matter inho- 
mogeneities outside the local NBG region. 

By tracing the mass distribution further out to a dis- 
tance of 100/i _1 Mpc, invoking the matter distribution 
in the complete PSCz mimicking mock samples, we will 
then evaluate the extent to which matter inhomogeneities 
within a 100/i _1 Mpc scale are able to account for the mo- 
tions within the local 30/i _1 Mpc region. From this we can 
infer in how far the external influence over the local region 
can be ascribed to matter fluctuations situated between a 
radius of 30/i _1 Mpc and 100/i _1 Mpc. 

In this study we also have to take into account the 
fact that a single 30/i _1 Mpc region cannot be considered 
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representative for the whole Universe, and generic con- 
clusions on the basis of the kinematics within a single 
30/i -1 Mpc volume cannot be drawn. This is also true for 
the the NBG mock samples in this work, even though they 
were selected according to some strict criteria (see section 
5.2). Analysis and conclusions will therefore be based on 
a straightforward average over the 10 different 30/i _1 Mpc 
mock samples which were constructed for each cosmolog- 
ical scenario. The dispersion in the extracted parameter 
values will provide a reasonable estimate for their signifi- 
cance. 

5.3.3. Analysis of Reconstructions 

The basic product of the FAM reconstructions are veloc- 
ity maps, in essence a velocity vector at the location of 
each galaxy in the sample. Our analysis consists of three 
different but complementary tracks. The first and most 
straightforward one is the visual inspection of the result- 
ing velocity vector maps. It provides a direct impression of 
the extent to which a FAM reconstructed field reproduces 
the true velocities. Also, it will provide a direct impression 
of a spatial coherence in the differences between true and 
reconstructed field, which is an incisive way to uncover 
systematic contributions like e.g. a bulk flow component. 

The second examination is a strictly local analysis, 
a pure point-to-point comparison between the velocities 
predicted by the FAM reconstructions on the one hand 
and the "true" A-body velocity of the same object on 
the other hand. To some extent, the analysis by means of 
scatter plots is the most direct and objective quantitative 
comparison between two fields. Various velocity related 
quantities will be assessed in this fashion. Note that these 
localized comparisons cannot address the presence of spa- 
tial coherence in the cosmic flows (even though they may 
uncover systematic effects caused by external influences). 

Finally, the third track is targeted towards a factual de- 
scription of the spatial coherence within the velocity fields 
or, rather, in the residual fields between the "true" veloc- 
ities and the reconstructed velocities. Systematic trends 
in these residual fields are interpreted as manifestations 
of external forces. Of these we shall determine the first- 
order - bulk flow - and second order - velocity shear - 
components. 

6. FAM velocity vector maps 

For reasons of consistency and to achieve optimal trans- 
parency the illustrated velocity vector maps in the follow- 
ing discussion all concern the same mock sample of NBG 
calculations. For the illustration of the FAM30 (Fig. 7) 
and the FAM100 (Fig. 8) reconstructions we use one of 
the ACDM 30/i _1 Mpc NBG mock catalogs. It is the same 
galaxy sample that was shown in 3-D in Fig. 5 and in 
projection along the "x-y" plane in Fig. 6. 

The vector maps in Figure 7 and Fig. 8 depict the 
projections of the raw unsmoothed galaxy velocities, for 
galaxies within a central slice of 10/i _1 Mpc. The size of 



the arrows is proportional to the amplitude of the pecu- 
liar velocity component within this slice, each arrow start- 
ing at the location of the galaxy. Both figures consist of 
three successive rows. The velocity maps in the first row 
correspond to the "real" world of the A-body simulation. 
The second row depicts the velocity maps for the FAM 
reconstructions, the FAM30 reconstruction in Fig. 7 and 
the FAM100 reconstruction in Fig. 8. The last row shows 
the resulting residual velocity vector fields, 

Vros = V F AM - V N body , (20) 

the vector difference between the A-body velocities and 
the corresponding FAM velocity reconstructions, [A-body 
- FAM 30 ] and [A-body - FAMioo]- 

Each row has three panels, containing the vector maps 
in the three mutually perpendicular "central" slices. Each 
plane is identified by means of the index combination "x- 
y", "x-z" or "y-z" (top figure), the index pair identifying 
the horizontal and vertical axis along which the panel is 
seen. Imagining these three planes passing through the 
centre of the 30/i _1 Mpc NBG volume provides a spatial 
impression of the full 3-D velocity field. Note that here the 
choice of Cartesian coordinate system does not have any 
special significance, arbitrarily set by the axes of the to- 
tal 345.6/i~ 1 Mpc simulation box (the "fundamental" box) 
from which the mock catalogs were distilled. This is unlike 
vector maps (e.g. Fig. 13) in some later sections. 

6.1. N -body sample: the "observed" velocities 

The velocity vectors in the top row vector maps depict the 
"real" A-body velocities of the "galaxies" located within 
the three "central" slices (the same for Fig. 7 and Fig. 8). 
The galaxy distribution is characterized by a few dense, 
massive and virialized clumps, visible as high concentra- 
tions of large and randomly directed velocity vectors. The 
trucly massive concentration visible in the lower left of 
the x — y panel is part of a superstructure extending be- 
yond the boundaries of the NBG region. It represents a 
major and dominant source for the motions in this area. 
This may be appreciated from the observed velocity flow 
towards this clump and the overall distortion of the flow 
in its vicinity. The large configuration visible in the "y- 
z" slice contains several dense compact regions embedded 
in a ridge-like structure running curvedly from the lower 
righthand corner to a location slightly left from the cen- 
tre. At least partially related to this mass concentration 
in and around the ridge is the bulk flow along the right- 
to-left direction. 

Overall, the "x-y" and "y-z" vector maps indicate the 
presence of a dominant coherent "bulk flow" pattern which 
can be traced throughout the whole NBG volume. By co- 
incidence, the orientation of the coordinate axes is such 
that the direction of the bulk flow is almost perfectly 
aligned along the "y"-axis: for this particular mock sample 
the "y"-axis does represent a physically significant direc- 
tion defined by the streaming pattern itself. The bulk flow 
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Fig. 7. 2D projected, unsmoothed velocities at each particles position in one of the NBG-ACDM mock catalogs. Each row 
contains three frames, corresponding to the central y — z, x — z and x — y plane through the point sample. Bottom row: 
projected iV-body velocities. Central row: projected FAM30 velocities. Bottom row: residual iV-body-FAM3o velocities. 



seems to be directed towards some (fictitious) point out- 
side the local 30/i~ 1 Mpc region. 

A dominant and conspicuous coherent flow pattern 
also characterizes the x — z velocity vector map. While 
the flow in the two other planes seems to be almost exclu- 
sively dominated by a bulk flow, here the pattern has a 
more complex geometry, readily recognizable as a typical 
"velocity shear" pattern. The specific shearing motion in 



this plane consists of a compressional component along the 
top lefthand to lower righthand direction, in combination 
with a dilational stretch along the perpendicular direction 
from the lower lefthand towards upper righthand corner. 



22 



E. Romano-Diaz et al.: External Tidal Impact on Local Supercluster Dynamics 



J. 




. .■ . ir . 

1 




j . 






- ft 




■ . ^ ■ f i 



■i ■ . 

■7 



■ 



■ " " " 

ir-L,.. 

Fig. 8. Same as Fig. 7, for the PSCz based FAMioo reconstructions. Shown are the 2D projected, unsmoothed velocities at each 
particles position in the central NBG mimicking volume of the full PSCz-ACDM mock catalogs, the same region as depicted in 
Fig. 7. Each row contains three frames, corresponding to the central y — z, x — z and x — y plane through the point sample. 
Bottom row: projected TV-body velocities. Central row: projected FAMioo velocities. Bottom row: residual TV-body-FAMioo 
velocities. 



6.2. NBG samples: FAM^q velocity vector maps 

The role of the local cosmic matter distribution on the 
motions in the local Universe is assessed on the based of 
the "FAM30 velocities". They are the peculiar velocities 
computed by FAM on the basis of the local matter dis- 



tribution, supposedly reflected by the galaxies within the 
NBG catalogs. The corresponding reconstructed velocities 
are shown in the second panel row of Figure 7. With their 
final position as boundary condition, each velocity vector 
is located at the same galaxy position as in the TV-body 
maps (top row). Note that the vector maps in Figure 7 



E. Romano-Diaz et al.: External Tidal Impact on Local Supercluster Dynamics 



23 



and Figure 8, and also the later ones in Fig 12, Fig 13 and 
Fig 14, show the pure unsmoothed velocity vectors (and 
do not "correct" for the virialized regions). 

The FAM 30 velocity maps are distinctly different from 
the corresponding TV-body velocity maps (top row): a co- 
herent flow pattern is almost entirely absent. The FAM30 
reconstructions obviously did not recover the strong bulk 
flow observed in the iV-body velocity maps, nor the strik- 
ing shear pattern in the x — z plane. Because the FAM30 
velocity field reconstructions solely relate to the matter 
distribution within the inner 30/i _1 Mpc NBG region, this 
indicates that the major share of coherent bulk flow and 
the velocity shear are due to the matter distribution out- 
side the central 30/i _1 Mpc. This is most readily apparent 
in the velocity residual maps [A-body - FAM 30 ] , the differ- 
ence between the TV-body and the FAM 30 velocity vector 
fields (bottom row of Fig. 7). In the residual field [iV-body 
- FAM30] we recognize the same characteristic flow pat- 
terns, strong spatial correlation, long-range coherence and 
overall morphology as in the full A-body velocity field. 
This represents convincing evidence for the external ori- 
gin of the large-scale "bulk" and "shear" component in 
the local velocity flow. 

Prominently visible in the residual velocity field is the 
strong bulk flow along the "y"-axis. Overall, the spatial 
pattern of the residual bulk flow appears to reproduce 
that of the A-body flow field. However, some minor yet 
significant differences between the residual and the full 
A-body bulk flow can be discerned. The amplitude of the 
corresponding velocities in the residual map is somewhat 
smaller than the equivalent A-body velocities: apparently 
part of the bulk flow is induced by the local NBG matter 
distribution. This docs not seem to be true for the ve- 
locity shear: the shear patterns in the "x-z" plane of the 
residual and iV-body velocity fields are almost identical 
(except for the virialized motions in high-density clumps). 
Apparently, the velocity shear component is almost exclu- 
sively due to external matter distribution. As a locally 
flattened matter configuration would induce an internal 
shear flow, this appears to imply a local matter distribu- 
tion whose geometry is hardly flattened or elongated. 

Closer inspection of the FAM30 velocity field provides 
a more detailed view of the small-scale flow pattern men- 
tioned above. In the "x-y" plane the large-scale (A-body) 
bulk flow has virtually completely disappeared. Instead, 
the dominant motion in the "x-y" plane is a streaming 
flow towards a prominent matter concentration within this 
region (lower left). On the other hand, in the "y-z" slice 
a trace of the TV-body bulk flow along the "y" axis re- 
mains, be it that the corresponding velocities have con- 
siderably smaller amplitudes than their A-body counter- 
parts. These local motions appear to be effected by the 
matter located along the lower ridge, supporting the im- 
pression that this feature is a local extension or outlier of 
the large-scale matter configurations responsible for the 
full bulk flow. Examination of the panels in Fig. 5 and 
Fig 6 indeed seems to suggest that the density ridge in 
the lower half of the "y-z" plane is indeed connected to 



structures just outside the NBG volume, while this per- 
haps may be true for the massive matter clump in the 
"x-y" plane too. This may not come as a surprise: the lo- 
cal matter distribution will to some extent be correlated 
with the external matter configuration so that the locally 
induced bulk flow is expected to reflect at least partially 
the full TV-body bulk flow. 

In summary, the inability of the FAM30 reconstruction 
to recover the large-scale bulk flow and velocity shear is 
a consequence of the fact that they are a result of the 
action of the mass distribution on scales larger than the 
internal 30/i _1 Mpc size region while the FAM30 velocities 
are entirely and self-consistcntly determined by the mass 
distribution within this interior region. The residual [A- 
body - FAM30] maps, which are a model for the possible 
findings of a real- world observational campaign, provide 
the most elucidating illustration of their "external" ori- 
gin. Even though they do provide convincing evidence for 
their external nature, they do not provide sufficient infor- 
mation to infer the identity and nature of the main source 
of the flow patterns. In principle, however, we may de- 
duce a substantial amount of information on the basis of 
a careful quantitative analysis: the work by Lilje, Yahil 
& Jones (1986) still sets a prime example. To this end, 
we will investigate the external matter distribution in the 
PSCz 100/i _1 Mpc sized regions. 

As a final note, we point to the rather artificial nature 
of velocity vectors in the vicinity of the massive clump in 
the "x-y" slice as indicative for the self-consistent nature 
of the FAM reconstructions. Its location near the edge 
of the NBG volume even appears to have generated the 
rather contrived infall motions along the rim of the NBG 
sphere. 

6.3. PSCz samples: FAMiqq velocity vector maps 

The contribution by the relatively nearby external mat- 
ter agglomerations, within a distance of 100/i _1 Mpc, to 
the motions in the local Universe is investigated on the 
basis of the "FAM100 velocities" . FAM produces these pe- 
culiar galaxy velocities on the basis of the galaxy sample 
in the full mock PSCz galaxy sample, extending out to 
lOO/i Mpc around the center of our local region. The 
corresponding reconstructed velocities are shown in the 
second panel row of Figure 8. It is the analogy for the 
"FAM100 velocities" of Figure 7, and concerns the same 
30/i _1 Mpc central region (the NBG region is the central 
subregion of the PSCz mimicking catalog). The FAM 30 
maps showed the dominant influence of externally induced 
forces on the motions in the local 30/i~ 1 Mpc NBG region: 
on the basis of the FAM100 maps we seek to assess whether 
the major share of the responsible external matter agglom- 
erations may be identified within the realm of a PSCz like 
volume. 

Comparison of the first and the second row of panels 
in Fig. 8 shows the large degree of similarity between the 
FAM100 velocities (panels 2 nd row) and the TV-body veloc- 
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ities (panels top row). Unlike the FAM 30 maps in Fig. 7 
we find that the FAMioo maps successfully reproduce most 
of the large-scale behaviour and most of the finer details 
of the A-body velocity field. The degree of similarity is 
particularly evident in the corresponding residual velocity 
field [A-body - FAM 30 ] (bottom row panels). With the ex- 
ception of the high-density virialized regions the residual 
velocities are very small and mostly randomly oriented: no 
significant spatial correlations and spatial coherence can 
be detected. 

The detailed similarity between the A-body and the 
FAMioo maps shows that it is sufficient to take account 
of the mass distribution out to 100/i _1 Mpc for explaining, 
in considerable detail, the velocity flows in the local NBG 
volume. Moreover, the detailed rendering of the velocity 
field by FAM is a convincing demonstration of the capac- 
ity of the FAM technique to accurately describe the dy- 
namics implied by the observed local galaxy distribution. 
The quantitative comparisons in the following sections will 
provide ample support to this claim. 

Of course, the above conclusion is partially related to 
the realizations of the cosmological scenarios we have stud- 
ied. The behaviour of the power spectrum P{k) on large 
scales will considerably influence the generality of our find- 
ings. A power spectrum with more power on large scales 
would modify our findings: potentially it may be so that 
we need a representation of the matter distribution out 
to larger radii than 100ft-^ 1 Mpc. In this respect it is im- 
portant to note that the used A-body velocity fields do 
not have any contributions from wavelengths larger than 
« 175h- 1 Mpc (both for tCDM as well as ACDM simula- 
tions, see Fig. 2). This merely for the technical reason of 
the simulation box imposing an upper limit to the scale 
on which we can represent P(k). The extent to which this 
may influence our conclusions may be readily appreciated 
from Figure 2 (right column, top and bottom panel: com- 
pare solid lines with dashed ones). The velocity field per- 
turbations of tCDM and ACDM carry out considerably 
further than the fundamental scale of the simulation box, 
in particularly affecting the resulting bulk flows. 



7. Point-to-point comparison 

Scatter diagrams are used to assess the point-to-point 
comparisons between quantitative aspects of the "real" 
galaxy velocities in the original A-body samples and the 
computed velocities in the FAM 30 and FAMioo reconstruc- 
tions. This analysis is meant to be a direct, in principal 
local, assessment of systematic trends in the velocity flows 
in volume of the NBG sample. The comparisons involve a 
component of the "true" A-body velocity (abscissa) versus 
the equivalent quantity for either the FAM 30 or FAMioo 
velocities, or of the corresponding residuals (ordinate). 
Since the objects had been artificially added to the Zone of 
Avoidance, any particles inside this region (see sec. 5.2.2) 
are excluded from these diagrams. 



7.1. Systematics 

If we neglect the small-scale sources in the deficiencies of 
FAM reconstructions, the differences between FAM 30 and 
FAMioo scatter plots are mainly to be ascribed to the 
corresponding differences in the external gravitational in- 
fluence acting over the two corresponding sample volumes. 
In the external gravitational influence the the correspond- 
ing leading velocity terms are the bulk flow Vbuik and the 
velocity shear Sij, 

3 

«Nbody,i « VFAMA + «bulk,i + Sij XX j + ... (21) 

j=l 

In the above i, j denotes the Cartesian component indices. 
The vectors Xj represent the vector components along the 
Cartesian j direction of the spatial unity vector oriented 
along the object position vector x. Following this defini- 
tion, xx j is the j-component of the position vector x, with 
x the distance of the object. 

Systematic differences in FAM velocity- A-body veloc- 
ity scatter diagrams are therefore to be attributed to dif- 
ferences in bulk flow, shear and possibly higher order con- 
tributions. Because each of these large-scale phenomena 
will manifest themselves in distinctly different ways, we 
seek to identify them from the scatter diagrams. An hor- 
izontal offset in the scatter diagram would be the trade- 
mark for a bulk flow component. Velocity shear would 
manifest itself as a distinctly characteristic correlation be- 
tween residuals and velocities, although the prominence 
of this signal will be dictated by shear magnitude, con- 
figuration, and orientation with respect to the reference 
system (as is true for the bulk flow). In reality, the situ- 
ation will be more intricate. Subtle correlations between 
small-scale and large-scale contributions will bring about 
a change in the slope of the scatter diagram of FAM re- 
constructed velocity components against their full A-body 
values (see 7.3.2). 

7.2. Velocity Scatter Diagrams Analysis 

Scatter diagrams are presented in three successive figures. 
The depicted scatter diagrams all relate to a ACDM mock 
catalogue, and each of these point-to-point analyses re- 
lates to a different aspect of the velocity field reconstruc- 
tions. Figure 9 contains four different panels, of which 
each contains two scatter diagrams: FAM 30 versus A-body 
quantity (left) and the equivalent FAMioo versus A-body 
quantity (right). The diagrams in Fig. 10 focus on the cor- 
relations between these quantities and the scatter around 
regression relations. The figure addresses three velocity- 
related quantities, each taking one column of each 2 pan- 
els: the top one for the comparison of the FAM 3 o compo- 
nents with their A-body counterparts, and the same for 
the FAMioo components in the bottom frame. 

A straightforward comparison is that between the 
Cartesian velocity components Vi of the FAM recon- 
structed velocities and the A-body velocities (Fig 10, 
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Fig. 9. Point-to-point comparison (scatter plot) of four quantities related to FAM reconstructed velocities and the corresponding 
"real" object velocity in the underlying JV-body samples. The iV-body realization, and the mock sample, concern the ACDM 
scenario. For each quantity two panels are shown: the lefthand one is the one for FAM30 reconstruction, the righthand one for 
the FAM100 reconstruction. Top lefthand: the amplitude of the velocities. Top righthand: residual velocity amplitudes FAM- 
iV-body. Bottom lefthand: FAM velocity component projected along the A^-body velocity. Bottom righthand: angle 9 between 
FAM velocity and -/V-body velocity, in terms of fi = cos(#). 



righthand column). Complementary regressions involve 
coordinate system independent aspects of galaxy veloci- 
ties. These involve the velocity amplitude |ufam| (Fig. 9, 
top lefthand panel), the component of each FAM veloc- 
ity parallel to the corresponding iV-body velocity, v\\ (= 
v proj), and the additional perpendicular component v± 
(Fig 10, first column). Misalignments between the real 
iV-body velocity and the FAM velocity reconstructions 
should indicate in how far a reconstruction has been fail- 
ing to take into account all relevant gravitational forces 
along the path of a particle. Systematic misalignments re- 
veal themselves in the scatter diagram of the angle 9 be- 
tween the FAM velocity and the galaxies' iV-body velocity 
VNbody (in Fig. 9 we plot fi = cos(0), bottom righthand 
panel). In terms of the character and systematics of the 
underlying physics and dynamics the residual velocities, 
[A^-body - FAM], represent highly informative aspects in 
our analysis. They are assessed in (Fig. 9, top righthand 
panel) and Figure 11. 

The significance and strength of correlations between 
the A^-body and FAM velocity components in the scat- 
ter diagrams are analyzed by means of a linear regression 
and correlation analysis. To circumvent excessive pollu- 
tion of the computed parameters by the virialized mo- 



tions in high-density regions (see sect. 3.3), the galaxy 
velocity components in these regression analyses involve 
2/i _1 Mpc tophat filtered velocity fields. The resulting nu- 
merical values of the correlation parameters are listed in 
Table 3. Table 3 is organized in two separate sections, 
one for the regression analysis results of the ACDM mock 
samples (top section) rCDM samples (bottom section). 
For both the ACDM and the tCDM section we list the 
results for four velocity related quantities, each separately 
for the FAM30 and the FAM100 reconstructions. 

The presence of significant correlations between FAM 
reconstructions and their A^-body counterparts is evalu- 
ated on the basis of the nonparametric Spearman corre- 
lation coefficient Rs P ear ■ The linear regression parameter 
Rim quantifies the linearity of the relation. Prevailing in 
most situations, the linear regression parameters are used 
to characterize the relation between reconstructed and real 
A^-body velocities: the zero-point (offset) arj, the slope ai rg 
and the dispersion ai rg around the linear regression rela- 
tion. We assume equal errors in FAM and iV-body ve- 
locities, as both are affected by similar shot noise errors 
(while 2ft, _1 Mpc top hat smoothing significantly reduces 
the impact of virial motions on FAM velocity predictions) . 
In addition, we also list the rms scatter of the parameters, 
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Table 3. Results of the linear regressions in the tests with the mock catalogs. Col. 1: Cosmological model. Col. 2: Component 
Col. 3: label of the set. Col. 4: the average Spearman (non-parametric) correlation coefficient Col. 5: average linear correlation 
index of the best fitting line and its la scatter. Col. 6: average zero point of the best fitting line and its la scatter. Col. 7: 
average slope of the best fitting line and its la scatter. Col. 8: average dispersion around the fit and its la scatter 



Cosmology 


Component 


bet 


D 


r> 


a 

(km s- 1 ) 


&lrg 


airg 

(km s- 1 ) 




\v\ 


FAM30 

FAMioo 


0.54 ±0.11 
0.68 ±0.12 


0.56 ±0.11 
0.62 ±0.16 


-29.43 ± 84.28 
83.38 ± 156.39 


0.66 ±0.21 
0.85 ±0.26 


198 ± 73 
212 ± 89 


ACDM 


v \\ 


FAM30 
FAMioo 


0.43 ±0.14 
0.59 ±0.23 


0.37 ±0.20 
0.52 ±0.29 


-198.21 ±173.86 
-135.98 ±242.78 


0.65 ±0.30 
0.99 ±0.40 


251 ±114 
281 ± 158 




Vi 


FAM30 
FAMioo 


0.54 ±0.20 
0.71 ±0.27 


0.48 ±0.29 
0.61 ±0.31 


19.35 ± 140.90 
5.96 ± 78.94 


0.66 ±0.25 
0.86 ±0.23 


240 ± 130 
233 ± 107 




Vres,i 


FAM30 
FAMioo 


0.65 ±0.11 
0.40 ±0.14 


0.66 ±0.11 
0.47 ±0.19 










\v\ 


FAM30 
FAMioo 


0.56 ±0.05 
0.66 ±0.10 


0.59 ±0.05 
0.59 ±0.13 


-54.10 ± 67.09 
-146.58 ± 140.40 


0.84 ±0.18 
1.34 ±0.26 


203 ± 57 
289 ± 80 


rCDM 


V \\ 


FAM30 
FAMioo 


0.52 ±0.09 
0.63 ±0.14 


0.48 ±0.15 
0.55 ±0.19 


-313.16 ± 134.10 
-430.59 ±302.29 


1.04 ±0.24 
1.60 ±0.46 


272 ± 100 
360 ± 171 




Vi 


FAM30 
FAMioo 


0.68 ±0.13 
0.79 ±0.17 


0.63 ±0.15 
0.69 ±0.13 


-4.75 ± 107.95 
27.60 ± 70.70 


0.80 ±0.19 
1.12 ±0.11 


221 ± 81 
255 ± 87 




V r es,i 


FAM30 
FAMioo 


0.55 ±0.15 
0.22 ±0.13 


0.57 ±0.13 
0.27 ±0.14 









estimated on the basis of the results for the 10 different 
mock catalogs (for each of the four different configura- 
tions). 

7.3. Inventory 

7.3.1. Velocity Amplitude 

In the top lcfthand frame of Fig. 9 the FAM30 and 
FAMioo velocity amplitudes are compared with their TV- 
body counterpart «Nbody|- The FAM30 diagram differs 
considerably from the FAMioo diagram: the FAM30 veloc- 
ities are systematically smaller than their FAMioo coun- 
terparts. Also, while the latter have a strong one-to-one 
correlation to the TV-body velocities, the FAM 30 diagram 
shows a systematic offset with respect to this relation (the 
solid line) and a somewhat larger scatter. 

While the FAM 10 o diagram tapers out to higher ve- 
locities and even shows a few points with |vfam| > 
2000 km s _1 , there is a firm ceiling of «fam ~ 
1300 km s _1 for the FAM30 velocities. It is a direct re- 
flection of the FAM30 reconstructions missing out on the 
gravitational force contributions by the external mass dis- 



tribution. The asymmetric nature of the scatter in both 
diagrams is due to particles in high density regions. 

To a good approximation, the correlation between the 
FAMioo velocity amplitudes and «Nbody| is that of a linear 
identity relation: the solid line, |t>Nbody| = I^famI, forms 
a good fit to the scatter diagram (see Table 3: ai rg — > 1). 
The significantly higher value of Spearman's correlation 



coefficient (Table 3: R,s 



pear 



0.68 vs. R s 



pear 



0.54 for 

FAM30) indicates and confirms the visual impression of 
Fig. 8) of the tight correspondence between the FAMioo 
and TV-body vector velocity fields. The FAM30 results 
stand in marked contrast: the majority of the FAM30 ve- 
locities have a systematically lower amplitude than their 
TV-body counterparts. It results in a relation with a sig- 
nificantly shallower slope than that of the identity rela- 
tion |«Nbody| = |«fam| (also see table 3): objects with a 
higher velocity have a larger discrepancy. The contribu- 
tion by the missing large-scale velocity component vi ss to 
the amplitude of the FAM velocity includes a cross-term 
(— VNbody • viss), a term dependent on the velocity VNbody 
of the galaxy. Most of the missing large-scale velocity com- 
ponent viss is due to the absence of a bulk flow term in 
the FAM30 reconstructions. Subtle and/or higher order 
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Fig. 10. Point-to-point comparison (scatter plot) of FAM velocity decompositions wrt. the corresponding TV-body velocity of the 
involved ACDM mock sample. Both the parallel projection vu (centre frames) as well as the perpendicular projection v± along 
the iV-body velocity direction are analyzed. Top row: FAM30 reconstruction against A-body velocity equivalent. Bottom row: 
FAM100 reconstruction against iV-body velocity equivalent. Lefthand column: Vu. Central column: v±. Righthand column: The 
x component of the velocity, v x . The solid line in the righthand column frames indicates the identity line v x ,fam = v x,N-body 
The dashed line indicates the best linear fit. Evidently, the FAM30 deviate systematically from the identity line. The point 
scatter density is characterized by the contours, indicating iso-density surfaces. See Table 3. 



external gravitational effects play an additional role: the 
velocity vector diagrams did already reveal that the pres- 
ence of shear should be one of the main contributors (cf. 
eq. 21). 

When comparing the ACDM FAM velocity amplitudes 
with those of the rCDM reconstructions it is evident that 
in the case of the FAM30 reconstructions the latter adhere 
considerably better to the corresponding iV-body values. 
The linear fitting slope ai rg (see table 3) is considerably 
closer to unity for the rCDM samples than for the ACDM 
samples. Over a 30/i _1 Mpc volume the external density in- 
homogeneities in the ACDM cosmology will induce consid- 
erably higher bulk flows than the more moderate rCDM 
perturbations, which is entirely in line with the theoreti- 
cal expectation (fig 2). In the case of the FAM100 recon- 
structions the qualitative differences are far less promi- 
nent. On the scale of ~ 100/i~ 1 Mpc the mass distribution 
in both the rCDM and ACDM simulation volumes have 



converged to homogeneity and no major bulk flows are to 
be expected. 

7.3.2. Velocity Decompositions 

In the two righthand frames of Fig. 10 we show the scat- 
ter diagrams for the x-component of the FAM30 (top) and 
FAM100 velocities (bottom). Although with a significant 
level of scatter, the FAM100 diagram can be fitted quite 
well by a straight line with a slope close to unity (linear 
regression line: dashed, unity line: solid). That the equiva- 
lent FAM30 diagram may also be fitted by a straight line, 
be it with a slope significantly smaller than unity is not 
entirely straightforward. It stems from an intricate inter- 
play between the small scale velocity field and its larger 
scale contributions, which in most circumstances are not 
uncorrelated (Berlind, Narayanan & Weinberg 2000). 
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Table 3 lists the linear regression parameters. 
Although the average best fitting slopes for the FAMioo 
velocities are either larger (rCDM) or smaller (ACDM) 
than unity, the deviation from unity is considerably 
smaller than that for the FAM 30 velocities, well within the 
~ lcr uncertainty interval. In all, these regression results 
do adhere to the expected and noted trend of FAMioo 
velocities accounting for practically all contributions to 
the local velocity field and FAM30 velocities systemati- 
cally neglecting significant external contributions. Notice 
that the scatter around the regression lines for the FAM 30 
and FAMioo reconstructions is of comparable magnitude 
(as may be inferred from the superposed number density 
contours). 

While the choice of any Cartesian coordinate system 
is an arbitrary one we have also addressed the decomposi- 
tion of the particle velocities in one defined by the system 
itself. The FAM velocities are decomposed in a component 
projected along the corresponding A-body velocity, v pro j 
(or v\\) and the complementary perpendicular component, 
v±. The second column of figure 10 contains the scatter 
diagrams for the parallel component of the FAM 30 (top) 
and FAMioo (low) velocities. Qualitatively, the behaviour 
of both diagrams resembles that of the velocity amplitude 
scatter diagrams in fig. 9. A 1 — 1 relation between FAMioo 
velocities and A-body velocity amplitude represents a rea- 
sonable fit (solid line: slope ai rg ss 0.99). The FAM30 dia- 
gram not only appears to deviate strongly from such a 1 — 1 
relation, it may even fail to fit any linear relation. Also, 
none of the projected FAM30 velocity components appears 
to supersede a value of ~ 1200 km s _1 . Given the fact that 
the equivalent FAMioo component even surpasses values 
of « 2000 km s _1 , this confirms the systematic deficiency 
of the gravitational field in the FAM30 evaluations. 

From the scatter diagrams for the perpendicular FAM 
velocity components, v±, one can infer that almost all sys- 
tematic effects are confined to the parallel components 
vn. For both FAM30 and FAMioo the complementary 
perpendicular component lacks a systematic correlation 
with the A-body velocity. It mainly represents unrelated 
scatter, with a magnitude concentrated around values of 
~ 200 — 250 km s _1 . The only difference between the 
FAM30 and FAMioo reconstructions is that for the lat- 
ter v± involves considerably higher values, reflecting the 
higher amplitude of the FAMioo velocities. FAM30 veloc- 
ities, on the other hand, involve stronger misalignments 
(section 7.3.3). 

7.3.3. Velocity Alignments 

Misalignments between the reconstructed FAM velocities 
and the A-body velocity vectors are the result of a few 
effects. A major source is that of localized small-scale ef- 
fects. These are not expected to lead to systematic offsets: 
they will have a noisy character and reflect random mo- 
tions in highly nonlinear environments, in particular those 
of dense virialized regions. Because these have no preferred 



direction, they behave like randomly oriented "residual" 
velocities wrt. to the real A-body velocities of galaxies. Of 
an entirely different nature are misalignments stemming 
from the systematic neglect of the external gravitational 
forces. Because the resulting residual velocity vectors com- 
prise systematic components along one or a few preferred 
directions, a distinctly anisotropic distribution is the re- 
sult. This reflect itself as a systematic trend for total N- 
body galaxy velocities to be aligned along the residual 
velocity components. 

For both the FAM 30 and FAMioo reconstructions we 
computed the angles between the A-body velocity vec- 
tor and the FAM velocities. The lower righthand panel 
in Figure 9 confirms that the alignments of "residuals" 
and total velocity is indeed considerably stronger for the 
FAM30 reconstructions than the FAMioo ones. The fig- 
ure plots, for each galaxy in the sample, the misalignment 
angle 8 (or, rather, /i — cos(#)) versus the A-body veloc- 
ity magnitude «Nbody For the FAMioo velocities we see a 
near isotropic distribution of angles. With the exception of 
a minor concentration near perfect alignment, cos(#) = 1, 
the distribution is sweeping out nearly uniformly over the 
full range of = cos(#) = 1 — > —1. If at all there is a trend 
in velocity amplitude, it appears to be the weak tendency 
for large velocities to be better aligned. 

The above results reflect the observation that FAMioo 
residual velocities do mainly consist of small-scale random 
effects. The FAM 30 residuals form a telling contrast. They 
are heavily aligned along the full A-body velocities, with 
a very strong concentration near 9 = 0. Although occa- 
sionally there are serious misalignments, their occurrence 
diminishes rapidly towards large 8. When they occur it 
almost exclusively concerns small velocities, mostly cor- 
responding to serious misalignments between the locally 
induced velocity and the added external velocity compo- 
nent. 

7.3.4. Velocity Residuals 

The residuals accumulate all systematic physical effects as 
well as random artifacts. They are therefore an excellent 
source of information on the dynamical role of matter con- 
centrations in the various galaxy sample volumes. If there 
are large external contributions to the galaxies' velocity 
these will constitute a major part of the residuals. On 
the other hand, if most of those influences arc contained 
within the sample volume treated by FAM, the residuals 
may mainly reflect localized nonlincaritics and artifacts of 
the FAM method. Scatter diagrams involving the resid- 
ual velocities will indicate systematic trends and are well 
suited for elucidating the character and underlying dy- 
namics of external influences. 

Figure 11 elaborates on this observation. In two succes- 
sive rows, the top one for the FAM 30 reconstructions and 
the bottom one for the FAMioo reconstructions, it displays 
the residuals v res for each of the three Cartesian velocity 
components, v x , v y and v z . Each panel plots the veloc- 
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Fig. 11. point-to-point comparison (scatter plot) of the Cartesian components of the residual velocity VFAM~VNbod y components 
against corresponding iV-body velocity "VNbody From left to right: v x , v y and v z . Top row: FAM30 reconstructed velocity. Bottom 
row: FAM100 reconstructed velocity. 



ity component residual as a function of the corresponding 
iV-body velocity component. 

The mark of a bulk velocity is a constant offset of the 
scatter diagram, a translation of all FAM velocities by 
a constant term. This is indeed what is observed in the 
FAM30 v y scatter diagram: the vast majority of points 
is located beneath the v y = km s _1 Mpc _1 line. It is a 
telling confirmation of the impression yielded by the corre- 
sponding velocity vector fields in Fig. 7. The velocity vec- 
tor field revealed the presence of strong bulk flow oriented 
almost perfectly along the y-axis: clearly visible in the N- 
body velocity field, hardly present in the corresponding 
FAM30 velocity field, and representing a major component 
of the residual field [Af-body - FAM 30 ] . When turning to 
the equivalent FAM100 diagram, the indicative offset for 
a bulk flow has almost completely disappeared. This im- 
plies that the source(s) for the bulk flow should be found 
within the region between 30/i _1 Mpc and lOO/i^Mpc. 
The equivalent v x and v z FAM 30 residual scatter diagrams 
do confirm the visual impression of there hardly being a 
bulk flow contribution along the x- and y-directions. 

Additional systematic behaviour is readily apparent 
in Fig 11: the diagrams show an almost linear increase 
of residual velocity with A^-body velocity. Also, we find 
that the v x scatter is skewed towards negative v res values 
while the v z diagram is skewed towards positive v res val- 
ues. In the equivalent FAM100 scatter diagrams the linear 



increase of v res and the asymmetry in the v x and v z dia- 
gram has almost disappeared: the dense core of points has 
turned into a compact and nearly horizontal bar symmet- 
rically distributed around the v res = km s _1 line. To a 
large extent this is explained by the much smaller contri- 
bution of external tidal shear to the flows over 100/i _1 Mpc 
volumes (cf. Fig. 2). 

The mark for external shear is a near linear increase of 
residual velocities as a function of their A^-body (or mea- 
sured) velocity. Depending on the location x of a galaxy 
within the sample volume and with respect to the shear 
configuration its participation in a shear flow will involve 
a velocity component v s = ^SijXj. This may involve a 
negative or a positive contribution. With such shear con- 
tributions representing a non negligible component to the 
total velocity, its systematic contribution to a largely ran- 
dom local residual signal reshuffles the velocities such that 
on average the largest velocity involves the largest residual 
contribution. 

With prominent large-scale bulk and shear motions at 
large, the FAM100 residual scatter diagram has largely 
transformed into a featureless and purely random point 
distribution. The residuals mainly involve uncorrelated 
small-scale effects and are nearly independent of the am- 
plitude of the A^-body velocity. Some additional artifacts 
are seen upon closer inspection: the presence of diffuse 
"S" -shaped point clouds in both the FAM 30 and FAM100 
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residual diagrams, tapering off towards a steep tail at both 
the negative and positive side of the plots. 

The corresponding scatter diagrams for the velocity 
residual amplitudes |v rcs | = |vfam — VNbodyl represents a 
summary of the systematic trends (Fig. 10). The FAM30 
velocity residuals show a near linear increase as a func- 
tion of the A-body velocity, starting with an offset, in- 
dicative of the ingredients of bulk and shear flow in the 
residuals. The lack of any clear correlation between |v res | 
and |v Nbody | in the case of the FAM 10 o residuals confirms 
the absence of such systematic components. More clearly 
than in the case of the individual Cartesian components, 
the presence of local nonlinear motions may be discerned 
from the extensive surrounding clouds of outliers. 

7.4. Power Spectrum Dependence 

The contrast between FAM 30 and FAM100 scatter dia- 
grams is more pronounced in the case of the ACDM mock 
catalogs than in those assembled for the rCDM universes. 
This clearly reflects the fact that within the ACDM sce- 
nario cosmic structure is characterized by a larger co- 
herence scale. It implies the presence of larger and more 
coherent structures whose size exceeds 30ft- _1 Mpc. Their 
combined gravitation impact will yield a stronger system- 
atic impact in the velocity- velocity comparisons. On the 
other hand, the dispersions listed in table 3 also show that 
it would hardly be possible to infer information on the 
cosmological scenario on the basis of one individual real- 
ization. The large dispersion around the average slopes, in 
particular in the case of the ACDM Universe, show that 
the magnitude of the external dynamical effects may vary 
appreciably as a function of the location of the (mock) 
NBC sample within the simulation box. Local measure- 
ments will therefore be unable to separate cosmological 
effects from those stemming from local variations. 

7.5. Nonlinearities 

The point-to-point diagrams discussed above all contain a 
substantial level of scatter around the inferred regression 
relations. With a few exceptions the scatter of velocity 
quantities is in the order of ~ 200 — 250 km s _1 , for both 
the ACDM as well as the rCDM FAM reconstructions. 
The main source for this scatter are the virial motions in 
the high density and mildly nonlinear environments. Also 
shot noise provides a substantial additional contribution. 
In the case of small filter radii, another source of scatter 
is formed by spurious very close pairs of points in the par- 
ent A-body catalog which for artificial reasons failed to 
collapse into a single object (Branchini, Eldar & Nusser 
2002). Scatter may also be due to higher order multipole 
components in the external gravity field. An inspection of 
the particle configurations and the velocity vector maps 
does unmistakably show significant systematic variations 
on top of dipolar and quadrupolar components. However, 
tests restricting the analysis to points in the central re- 



gions of the sample produced no substantial decrease in 
level of scatter. This seems to argue for a minor role of 
such contributions. 

8. Bulk Flow and Tidal Shear: 

Velocity Flow Multipole Components 

In the previous sections we have found that in order to 
obtain a good representation of the local cosmic veloc- 
ity field it is necessary to take into account the external 
gravitational influence. This was accomplished through 
the incorporation of the fully detailed external mass dis- 
tribution contained in the (flux-limited) galaxy catalogs. 
This involved the galaxy distribution out to distances of 
100/i _1 Mpc. The reconstructions showed that modelling 
of velocity fields by FAM with the inclusion of matter 
concentrations on such large scales is indeed rewarding. 

In nearly all situations where the local volume Vint is 
suitably large, the small-scale details of the external mass 
configuration are rather irrelevant for constructing an ap- 
propriate model of the flows in the local Universe. An ap- 
propriate approximate expression for the the gravitational 
potential <fr ext (r) inside the internal volume V int due to the 
surrounding external matter distribution follows from its 
expansion in multipole contributions. Assuming a spheri- 
cal local volume with radius Ri nt , the potential & ext may 
be written in terms of a multipole expansion of spherical 
harmonics Yi m (0, </>) (see e.g. Jackson 1975) 

(22) 

00 m—l . « 

= "E E ^QimY lm (0,4>)r l , 

1=0 m=-l 

in which the multipole moments Qi m relate to the external 
density field p(x') as 

/•OO 

Q lm = / p(xV- 3 ^n(0',^x'. (23) 

J r 

Most contributions to the external gravity g, ext will be 
confined to these dipole and quadrupole components, 
induced by the corresponding large-scale constellations 
in which the surrounding matter concentrations have 
grouped themselves. Here we will assess the approximation 
in which the potential expansion (eqn. 23) is restricted to 
the monopole term I = 0, the dipole term I = 1 and the 
quadrupole term 1 = 2, 

*ext(r) » *o(r) + *i(r) + * 2 (r) (24) 
with 

m—l . , y 

*»( r ) = - E ^QimYue^)^. (25) 

m— — I 

To explore the nature of the external component in the 
total gravitational field in the Local Universe we proceed 
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by probing it through the resulting peculiar velocity field. 
The amplitude of higher order terms may be assumed to 
be so small that one cannot expect to deduce any signifi- 
cant value, given the sizeable errors in the available galaxy 
peculiar velocity datasets, we may expect this to be a rea- 
sonable approximation. We investigate the velocity field 
by decomposing the residual velocity field - i.e. the com- 
ponent in the velocity field which could not be accounted 
for in the FAM30 reconstruction and supposedly induced 
by external influences - into its multipole components. 
Once we have determined the bulk flow component and 
shear tensor components in the tidal velocity field, we will 
assess whether we can indeed relate this to the external 
gravitational ( "tidal" ) influence within the local Universe. 

Restricting the description of the external gravita- 
tional influence to the first few orders of its multipole 
expansion has several advantages. The large-scale exter- 
nal dipole and quadrupole gravity perturbations retain a 
largely linear character, simplifying the velocity field anal- 
ysis and thus retaining the direct linear relation between 
gravity and velocity field. Also, by discarding its small- 
scale fluctuating contributions a physically more transpar- 
ent image of the velocity field is obtained. This allows a 
straightforward relation and translation towards the cor- 
responding large-scale pattern of the surrounding mass 
distribution. A final practical issue of some importance 
is the fact that the dipole and quadrupole characteriza- 
tion is particularly suited for an implementation in FAM. 
Restricting the external force field to these moments al- 
leviates the need to take into account a large sample of 
external galaxies. Not only is the latter computationally 
expensive, in practice it is even not always feasible. 

8.1. Velocity Field Multipole Decomposition 

In the multipole analysis we restrict ourselves to the ex- 
ternally induced velocity components, v ex t, which in the 
following we frequently designate by the term "tidal" 1 . 
For each object, the "tidal" velocity vector is determined 
by subtracting the internally induced velocity field, v int , 
from the object's full velocity. The latter is usually the TV- 
body velocity of the mock galaxy, although we will assess 
the possibility of using the FAM100 velocity as a reason- 
able alternative. The internal velocity v int is deduced by 
evaluating, through our FAM computations, the impact of 
the internal matter distribution within the internal cata- 
log volume Vint- The resulting (residual) peculiar veloc- 
ity vector v ext field may then be expressed in terms of a 
Taylor series description as function of spatial position x. 

For the practical implementation, we follow the gen- 
eral scheme described by Kaiser (1991). The velocity field 
Taylor expansion is truncated at the quadratic term and 
is restricted to the dipole and quadrupole moments (and 

1 in the following we regularly use the word "tidal" to shortly 
indicate the externally induced component of a gravity or ve- 
locity field. As it includes a dipolar contribution, strictly speak- 
ing this is not an appropriate term. Also see sect. 3.6.2. 
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a minor monopole term). The tidal velocity field v ex t, is 
then modeled by the the first two components, a bulk flow 
vector, Ui, and a quadratic shear tensor contribution, , 

v ex t,i = Ui + SijXXj, where i,j = {1, 2, 3} , (26) 

in which i,j denotes the Cartesian component indices. As 
in eqn. 21, the vectors x j represent the vector components 
along the Cartesian j direction of the spatial unity vector 
oriented along the object position vector x. Using these 
notations, we can easily reconfigure eqn. 26 and express 
the i-component of the velocity of object n into a product 
of the vectors F n j and Vn, 

4 

v n ,i =^2F n j(x)V Il 
7=1 

in which the data 4-vector Fj and the velocity field com- 
ponent 4-vector Vn are defined as 

F nJ = {l,xxi,xx 2 ,xx 3 } 
Vn = {ui,sn,Si2,s a } . (27) 

The "dipolar" bulk flow components Hi and "quadrupo- 
lar" velocity shear components can then be obtained 
by solving for the vectors Vn on the basis of a fitting 
analysis (to be precise, s%j also includes a minor residual 
"monopole" expansion/contraction term). We accomplish 
this by computing for each Cartesian component i the val- 
ues for the multipole elements Ui and §ij which minimize 

x 2 

N oije 2 

X 2 = ( v n,e X t,i-J2F n j(x)Vn) , (28) 

n=l I 

to be evaluated on the basis of the data sample of N a bj 
objects at locations x„ and with inferred "external" ve- 
locities v„ jea;t . The bulk flow and velocity shear in the 
externally induced velocity component v ext , along with a 
residual expansion term, 

Vext,i ~ Vexp.i ~t~ Vbulk,i ~\~ ^shear.i (^9) 

will follow directly from the inferred values of the 4-vector 
components Vn: 

Vexp,i g -^~ I ^*(^) XXi 

3 

Vbulk,i = &i ; v shear,i — ^ ' ^^jXXj (30) 

i=i 

in which Tr(s) is the trace of the tensor Sij and s'^ the 
traceless shear tensor 

Tr(s) = sn + §22 + S33 

~s'ij = ~ Si] -\tt(§)5 13 . (31) 
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Field 

Configuration 



Cosmology v buik 

(km s" 1 ) 



si 



S2 



S3 



ACDM 313 ±65 13.3 ±3.8 8.8 ±3.6 0.6 ±2.5 -9.4 ±2.3 (kms^Mpc^ 1 ) 

rCDM 216 ±57 13.6 ±4.6 8.9 ±3.1 1.1 ±2.0 -10.0 ±3.7 

ACDM 399 ± 115 264 ± 107 19 ± 76 -283 ± 70 (km s" 1 ) 

rCDM 409 ± 137 267 ± 92 33 ± 60 -300 ±110 



Af-body - FAM30 



A^-body - FAM100 



ACDM 
rCDM 

ACDM 
rCDM 



81 ±51 
121 ±54 



3.7 ± 1.2 

4.8 ±1.1 

110 ±36 
145 ± 33 



2.6 ±0.9 -0.2 ±0.8 
3.4 ±1.0 0.0 ±0.8 



78 ± 27 
101 ± 30 



-6 ±25 
-1 ±23 



-2.4 ±1.0 (km s" 1 Mpc" 1 ) 
-3.3 ±0.8 

-72 ± 30 (km s" 1 ) 
-100 ± 23 



ACDM 317 ±89 14.0 ±4.6 9.4 ±4.0 0.5 ±2.9 -9.9 ±3.0 (kms^Mpc^ 1 ) 

rCDM 245 ±60 15.3 ±4.3 9.8 ±2.5 1.4 ±2.9 -11.2 ±3.8 

ACDM 421 ± 138 281 ± 120 15 ± 86 -296 ± 91 (km s _1 ) 

rCDM 458 ± 128 293 ± 75 42 ± 87 -335 ± 115 



FAMi 



FAM 3 



ACDM 81 ±51 2.4 ±1.3 -0.6 ±1.7 -0.2 ±1.4 0.9 ±1.3 (km s" 1 Mpc" 1 ) 

rCDM 121 ±54 2.9 ± 0.8 -0.6 ± 1.6 -0.5 ± 1.6 0.9 ±1.7 



A-body - FAM mp ; 



Table 4. Average results for the tidal bulk amplitude and shear eigenvalues. The errors represent the la scatter around the 
average value. Col. 1: Tidal field region where the velocity components were computed. The first two sections refer to the tidal 
field computed from the Af-body - FAM models. In the third section the external tidal influences are determined by the large 
region modeled on FAM100 instead of the A^— body samples. The fourth region refers to the residual field between N— body and 
FAM mp ; after the external contributions have been accounted for in the FAM modeling: FAM30 ± tidal bulk ± tidal shear (see 
text sec. 8.2 and sec. 8.5 for details). Col. 2: Cosmological model. Col: 3: bulk flow, iWifc- Col: 4: shear eigenvalue amplitude, 
s|. Col: 5: si eigenvalue (stretching). Col. 6: S2 eigenvalue (middle). Col. 7: S3 eigenvalue (compressional). 



8.2. Velocity Multipole Analysis: results 

The results of our analysis are summarized in Table 4. 
It lists the average quantities for the tidal bulk flow and 
shear components for the two cosmological scenarios dis- 
cussed in this work. The table has been organized in four 
(horizontal) sections. Each corresponds to another "differ- 
ential" velocity field, the difference between two differently 
processed velocity fields. 

For both the ACDM and the rCDM model each of the 
quoted values in Table 4 involve the average and standard 
deviation determined on the basis of ten different realiza- 
tions. This adds up to 8 configurations, two cosmologies 
per section. For each of the 8 configurations, in the third 
column the table lists the dipole component of the exter- 
nal velocity field, the bulk flow Vbuik- Subsequently, the 
velocity shear is specified in terms of the three eigenval- 
ues s\, S2 and S3 of the traceless shear tensor. This is 
preceded in the fourth column by the amplitude s of the 
shear. Note that shear is quoted in two units. First, units 



of km s _1 Mpc _1 , followed by the equivalent velocity dif- 
ferential in km s _1 over a volume of 30/i~ 1 Alpc radius. 
The intention of the latter is to offer a directly appreciable 
comparison between the relative importance of bulk flow 
and shear contributions. 

Each of the four sections specifies the values of the 
computed dipole and quadrupole moments of the velocity 
field of the corresponding sample. The first section relates 
to a multipole analysis of the differential velocity field be- 
tween the full A^-body velocity field and the FAM 30 ve- 
locity reconstructions of the inner 30/i _1 Mpc region, N- 
body - FAM30. The resulting residual velocity field has 
been generated by the mass distribution beyond a radius 
of 30/i _1 Mpc. On these linear scales the inferred dipole 
and quadrupole components of the velocity field may be 
directly related to the moments of the surrounding mass 
distribution. 

The second section of Table 4 does the same for the 
larger 100/i _1 Mpc region. The outcome of similar analyses 
are presented in the third and fourth section. The third 
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Fig. 12. Velocity Field Multipole Decomposition: the total JV-body velocity field, involving the same one as in Fig. 13 and 
Fig. 14, decomposed into its four different components. The coordinate system is that defined by the tidal shear tensor, see 
Fig. 13. At each row we depict the velocities in the (x — y), (x — z) and (y — z) plane. Top: the locally induced velocity, 
approximated by the FAM30 realization. Second row: the monopole component, a result of the slight local expansion due to its 
underdensity wrt. the global Universe. Third row: the bulk flow, of which most is concentrated in the (x — y) and (y — z) plane. 
Fourth row: the shear flow component. 



section repeats the analysis of the first section, except that 
the external tidal influences are determined on the basis 
of the difference between the FAM velocity reconstruc- 
tions within the large 100/i~ 1 Mpc region and the inner 
30/i _1 Mpc region. Earlier, in Section 6.3, we have found 
that the major share of the origin of the external tidal 
field is confined to this region and that it therefore may 
well be determined from the residuals between FAM100 
and FAM30. The comparison between the inferred multi- 
pole moments of the velocity differences between FAM 100- 
FAM30 in the third section and those in the first section 
are therefore expected to be rather similar, any systematic 
differences originating in tidal effects generated beyond a 
radius of 100/i _1 Mpc. The fourth section in Table 4 refers 
to the values of the residual tidal velocity field between N- 
body - FAM mp / , the FAM sample after having accounted 
for the missing external tidal contributions, FAM 30 + tidal 
bulk + tidal shear (see sec. 8.5). If indeed all significant 
contributions can be characterized by their dipolar and 
quadrupolar contributions, the multipole values in this 
section are expected to be negligible. 

8.3. Velocity Multipole Contributions: Maps 

For a direct visual appreciation of the various multi- 
pole contributions to the tidal velocity field we assess the 
"tidal" velocity field N-ho&y - FAM 30 , the velocity field 
generated by the mass distribution beyond a radius of 
30/i _1 Mpc, for one of the ACDM catalogs. The presented 
maps concern the same ACDM catalog as those presented 
in the maps of Figures 7. The map of the projection of this 
"tidal" velocity flow onto three central planes is shown in 
the top row of Figure 13. 

8.3.1. Dipolar component: bulk flow 

The externally generated velocity flow is dominated by its 
bulk flow component. This is in general true for both cos- 
mologies. The large impact of the bulk flow over the local 
30/i _1 Mpc volume can be immediately inferred from the 
values in the first section of Table 4, revealing contribu- 
tions in excess of 200 km s _1 . 

To facilitate visual appreciation of this observation we 
have have reoriented the reference system in Figure 13 
such that the x-axis is oriented along the bulk flow. While 
the original Cartesian system is an arbitrary one and thus 
lacks a physical context, the "bulk flow reference system" 
confines the inferred bulk flow u exclusively to the in- 
direction. As a result there are no bulk flow components 
in the corresponding y- and z-direction (note that within 



the y — z plane their direction is arbitrarily defined) . The 
pre-eminence of the bulk flow component can be immedi- 
ately seen in the x — y and x — z frames in the top row 
of Figure 13. Note that the same velocity maps, mostly so 
the y — z frame, reveal a clear shear pattern. 

8.3.2. Quadrupolar component: velocity shear 

Seeking to assess the quadrupolar term in the external 
velocity field we first remove the remaining expansion 
term from §ij. Diagonalization of the resulting traceless 
shear tensor s J ■ yields the shear eigenvalues and eigenvec- 
tors. The eigenvalues 81, 82 and S3 are indicative for the 
strength of the tidal force field induced by the surround- 
ing matter distribution, while the principal directions of 
this quadrupolar velocity perturbation field are indicated 
by the corresponding eigenvectors e Sj ,. 

The "shear ellipsoid", the quadratic surface defined 
by the shear tensor s'^ with principal axes aligned along 
the eigenvectors and with axis size set by the correspond- 
ing eigenvalue Sj, defines a natural reference system to 
assess the tidal shear flow field. The coordinate axes of 
this "shear reference frame" are identified with the or- 
thonormal basis defined by the (normalized) eigenvectors. 
The x-axis is chosen to be aligned along the major axis of 
the "shear ellipsoid", the direction defined by the largest 
(positive) eigenvalue s\ and directed along the strongest 
dilational (stretching) motion incited by the external tidal 
field. Likewise the z-axis is chosen to coincide with the low- 
est (negative) eigenvalue S3, aligned along the strongest 
"compressional" component of the tidal velocity flow. This 
leaves the y-axis as the one coinciding with the interme- 
diate eigenvalue 82- 

The imprint of the shearing motions can be discerned 
within the y — z plane and, most prominently, along the 
"x-z" projection of the "bulk flow reference system". After 
subtraction of the bulk flow component, i.e. A-body- 
FAM3o-tWfc , the quadrupolar component of the exter- 
nally induced velocity flow represents its principal con- 
stituent (Fig. 13, lower row). This is confirmed by the val- 
ues quoted in Table 4 for the shear contribution. In partic- 
ular when stated in the velocity equivalent unit of km s^ 1 
these shear values suggest that the quadrupolar shear con- 
tributions are of a comparable magnitude to those of the 
bulk flow. The maps in the lower row of Figure 13 suggest 
that there are strong dilational and compressional motions 
within the y — z plane. By contrast, the shear motions in 
the x-direction appear to be uncommonly weak. Given 
the "bulk flow reference system", it implies that for this 
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Fig. 13. 2D projected peculiar (residual) velocities for the same mock catalog as in Fig. 7 and Fig. 8 in three perpendicular 
central planes of 10/i~ 1 Mpc width. The coordinate frame is rotated such that the bulk flow velocity is oriented along the in- 
direction. Within the plane perpendicular to the ir-axis, the y and z axes are chosen arbitrarily. First row: the residual velocity 
(ie. tidal velocity) N— body - FAM30. Second row: residual velocities after subtraction bulk flow component. The resulting 
residual field is clearly dominated by a shear pattern, most notably in the y — z plane. 



particular realization we see a bulk flow directed almost 
perpendicular to the shear flow motions. 

Figure 14 depicts the same ACDM mock sample as 
presented in Figs. 7, 8 & 13, here in the "shear reference 
frame". The top row shows the full externally induced flow 
field, iV-body - FAM30, in this reference system. The tidal 
shear flows are almost exclusively confined to the x — z 
plane. This is most evidently illustrated in the central row 
of frames showing the velocity field without its bulk flow 
component: hardly any systematic flow is noticeable in the 
y-direction of the intermediate shear eigenvalue. 

8.4. Multipole Scale Dependence 

When turning to the external influences over a large 
lOO/i Mpc region, we may conclude from the second sec- 
tion of the table that most of the external contributions 
are accounted for, both bulk flow and shear are at least 



a factor of 3-4 smaller than for the inner 30/i -1 Mpc re- 
gion. The third and fourth section show that the explicit 
contributions from the regions between 100/i Mpc and 
30/i _1 Mpc and those beyond 100/i _1 Mpc are indeed sig- 
nificantly different, those beyond 100/i _1 Mpc tending to- 
wards zero contributions and as far as the shear is con- 
cerned almost an order of magnitude smaller than the 
equivalent contributions by the 30 — 100/i _1 Mpc region. 

A similar graphical assessment involving the FAM100 
reconstructions emphasizes the minor significance of tidal 
contributions stemming from density fluctuations beyond 
a radius of 100ft _1 Mpc. No coherent velocity pattern can 
be recognized in the residual velocity field between full 
A^-body and FAM100 reconstruction. The comparison be- 
tween this residual velocity field with the velocity maps 
including the contributions of the inferred bulk flow and 
shear flow do hardly show any difference. In all cases the 
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Fig. 14. 2D projected peculiar velocities for the same mock catalog as in Fig. 7 and Fig. 8 in three perpendicular central planes 
of 10/i _1 Mpc width. The coordinate frame is defined by the three eigenvectors of the shear tensor. The x-axis is aligned along 
the direction of the largest (stretching) eigenvalue, the j/-axis is aligned along the middle eigenvalue and the z-axis along the 
smallest (compressional) eigenvalue. The first two rows are the same as Fig. 13, within the "shear reference system" . The third 
row depicts the final residual field after subtraction of the quadrupolar shear component. 

velocity fields are dominated by the same thermal mo- 8.5. Multipole Velocity Flow Model 
tions. 

Following our argument that the externally induced ve- 
locity flow within the inner 30/i _1 Mpc mainly consists of 
a bulk flow and shear contribution, we may expect that 



E. Romano-Diaz et al.: External Tidal Impact on Local Supercluster Dynamics 



37 




Fig. 15. point-to-point comparison between the three successive residual velocities and the corresponding mock catalog TV-body 
velocity. The three panels in each row correspond to the x, y and z velocity components. The coordinate system is the "shear 
reference system". Top row: the tidal velocity field TV-body- FAM. Middle row: residual velocity field after subtraction bulk 
flow. Bottom row: residual velocity field after subtraction of both tidal bulk and shear components. 



the effect of the external gravity field can be sufficiently 
accounted for by adding these components to a local ve- 
locity field model based on the mass distribution in and 
around the Local Superclusters. 

By separating the "internal" FAM velocity field from 
the "external" multipole contributions of the (monopole,) 
dipole and quadrupole components of the "tidal" veloc- 



ity field and adding the two, we obtain a total "FAM- 
multipole" model velocity Vf am _ mp i i, 

v fam—mpl,i ~ v FAM,i v exp,i ~T" v bulk,i "T" Vshear,i (32) 

A visual impression of the extent of the successive mul- 
tipole contributions may be obtained from Fig. 12. The 
vector plots of the four velocity contributions to Vf am - mp i 
(Eq. 32) are depicted in four successive rows, each within 
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Cosmology Rspear Rl in &lrg ^Irg 

(kins" 1 ) (kins" 1 ) 



ACDM 0.62 ± 0.17 0.58 ± 0.22 -4.03 ± 75.45 0.83 ± 0.20 246 ± 134 

rCDM 0.75 ± 0.11 0.71 ± 0.14 -9.31 ± 30.19 0.93 ± 0.13 219 ± 86 



Table 5. Average final results for the A— body velocities vs. the corrected FAM30 after adding to it the tidal Bulk and Shear 
contributions. The errors represent the la scatter around the average value. Col. 1: Cosmological model. Col. 2: non parametric 
(Spearman) correlation coefficient. Col. 3: linear correlation index. Col. 4: zero point of the best linear fit. Col. 5: slope of the 
best fitting line. Col. 6: dispersion around the fit. 



the mutually perpendicular three central slices (wrt. the 
shear reference system) . The top row concerns the FAM30 
velocity field reconstruction, followed successively by the 
expansion/contraction term (monopole), the bulk flow 
(dipole) and velocity shear (quadrupolar) . 

From fig. 13 and fig. 14 we conclude that the differ- 
ences between the "full" iV-body velocities and Vf a m-m P i, 
the total sum of the internal FAM30 and external dipole 
and quadrupole contributions, do not appear to show sys- 
tematic trends as it can be noticed from the residual bulk 
and shear components in section 4 of Table 4. Wherever 
there are large deviations, these are mainly confined to 
the high density virialized regions. 

8.6. Point-to- Point Comparison 

A quantitative quality assessment of the "FAM-multipole" 
model is offered by the point-to-point comparison between 
the full ./V-body velocity and its difference with respect 
to the successive modes of the "FAM-multipole" veloc- 
ity in Fig. 15. The v x , v y and v z of the various veloc- 
ity components refer to the "bulk flow reference system". 
The top row, plotting v^body vs. the residual A-body- 
FAM30 , reveals the expected systematic differences due to 
missing externally induced contributions. Given the fact 
that the bulk flow in this reference system is confined 
to the x-componcnt, we may note the uniform system- 
atic shift of the x residuals with respect to the zeropoint 
(vNbodyiVres) = (0,0) (top lcfthand frame). The subse- 
quent addition of the dipolar bulk flow contribution to 
FAM30 leads to a systematic downward uniform vertical 
shift of Af-body-FAAI 30 -Bulk (middle row Fig. 15): also 
the residuals in the x-direction now center on v res = 
(note that by virtue of the bulk flow the A^-body veloci- 
ties in the ^-direction are also skewed to values larger than 

VNbody = 0). 

The three point-to-point diagrams in the middle row 
of Fig. 15 show that even while the bulk flow is taken into 
account systematic motions remain in all three directions. 
The point-to-point comparisons still follow a strong corre- 
lation with respect to the the ./V-body velocities. It mainly 
involves the presence of the quadrupolar velocity shear 
component (in addition to a minor ingredient contributed 
by the monopole expansion/contraction term). This can 



be immediately inferred from the comparison between the 
diagrams in the central and lower row of Fig. 15: once 
the quadrupole component "Shear" has been added to the 
"FAM3o+Bulk" velocities the systematic effects seem to 
have largely vanished. What remains in the residuals is 
mainly random scatter, centered on the v res = km s _1 
line, with some exceptional outliers originating in the viri- 
alized regions. 

We have quantified the point-to-point comparisons by 
performing linear regressions similar to those presented 
in Section 7. Table 5 summarizes the results of this com- 
parison for all catalog samples for both cosmologies. In 
both cosmological models the slope of the best fitting line 
is consistent with unity at ~ lcr confidence level. As ex- 
pected, the scatter around the fit is similar to that of all 
previous analyses (see Table 3). Offsets around the zero- 
point are consistent with zero, although with a large dis- 
persion. The strength of the point-to-point correlations 
has increased considerably with respect to their FAM30 
counterpart (Table 3) and it is very similar to the FAM100 
case. 

8.7. Surrounding Matter Distribution: Tidal Source 

The surrounding external matter distribution is the source 
for the tidal velocity field which we inferred in the previ- 
ous sections. For various purposes we wish to relate the 
computed dipolar bulk flow and quadrupolar shear flow 
components to the surrounding matter distribution which 
induced them. 

The induced tidal velocities involve spatial scales rang- 
ing from 30/i _1 Afpc to 100/i~ 1 Mpc. Over this range the 
linear theory of gravitational instability holds to good ap- 
proximation. This translates into a direct linear relation- 
ship between induced velocity w ext and the cumulative 
external gravitational force g ex t, 

2/(n,A) 

V ex t(x,i) = ^ Eext{X-,t), (33) 

with /(Q, A) the linear velocity growth factor. This linear 
relationship also holds for every component of the velocity 
and gravity fields, and thus also for the individual dipolar 
and quadrupolar components of the externally induced 
velocity field. They are directly proportional to equivalent 
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Fig. 16. AitofF projections of the galaxy distribution of a NBG + PSCz ACDM catalog. The top-most panel shows the 
corresponding NBG distribution, d sur = [0, 30]/i~ 1 Mpc. The subsequent panels depict the external galaxy distribution enclosed 
by the shells defined by survey depth d sur = [30 — 55]fa _1 Mpc, d sur — [55 — 70]/i _1 Mpc d sur = [70 — 85]/i _1 Mpc and d sur — 
[85 - lOOj/i^Mpc. 



dipole and quadrupolar components of the gravity field: 
2/(0, A) 



Vfc ui fc(x,t) 



3#ft 



buZfc (x, t) 



Sij (x, i) — 3HQ ^ J ^ 



(34) 



with the (external) gravitational tidal shear tensor %j is 
defined as (see van de Weygaert & Bertschinger 1996), 



_ d 2 &tidal 

10 dxidx; 



(35) 



Notice that because of its external nature, the term 
^V 2 ^>udai $ij is always equal to zero. 

Ideally, we would like to infer the external tidal po- 
tential <&udai directly from the galaxy distribution in a 
sufficiently large surrounding region. This is specifically 
true for its dipolar and quadrupolar moments, with the 



intention to insert these terms directly into the expres- 
sion for the FAM potential (eqn. 12 and eqn. 10). The 
required externally induced bulk flow velocity and veloc- 
ity shear should be the result. The comparison of the FAM 
computed velocities for the local volume, in combination 
with the computed tidal velocities (Sec. 8.2), and the ob- 
served and measured velocities would then enable us to 
determine the amount of mass and average density in the 
local volume. 

To determine the gravitational influence of the sur- 
rounding matter distribution, we set out to assess the 
sky distribution of galaxies in the local NBG volume, 
out to tmbg = 30/i _1 Mpc, along with the external mass 
distribution in radial shells out to a distance rpscz < 
100/i~ 1 Mpc. A prominent dipolar matter configuration in 
the sky distribution will translate into a strong bulk grav- 
ity force. Similarly, quadrupolar anisotropies will translate 




Fig. 17. Sky distribution of galaxies in external shell of survey depth d sur — 30 — 100/i~ Mpc. This galaxy distribution should 
reflect the mass distribution inducing the local tidal velocity flow. By means of symbols we have indicated the track of the 
gravity dipole and quadrupole eigenvector directions (cf. eqn. 36 and eqn. 37) on the sky, by radially expanding outward the 
survey depth d sur in steps of lfo _1 Mpc, from d sur = 30ft -1 Mpc to d sur — 100/i _1 Mpc. For comparison the same symbol, but 
then enlarged, indicates the directions for the corresponding bulk velocity and velocity shear eigenvector directions. Diamond: 
dipole. Triangle: stretching component shear. Square: middle component shear. Star: compressional component shear. 



into an effective tidal shear force. In figure 16 we have plot- 
ted the galaxies in one of our ACDM mock catalogs in five 
successive distance shells. Aitoff projections of the angu- 
lar positions of the galaxies, as seen from the centre of the 
local NBG volume, provide an impression of the level of 
anisotropy in the mass distribution at successive radii. 

The first sky plot (top sphere) depicts the sky posi- 
tion of the galaxies in the local NBG-mimicking mock 
sample. It involves a highly flattened distribution, per- 
haps reminiscent of the Supergalactic Plane. The four 
subsequent shells correspond to successive cuts through 
the PSCz mimicking samples, at sampling depths d sur = 
[0-30], [30-55], [55-70] and [85-100] /i^Mpc. The first 
and direct observation is the diminishing sample density 
as a function of survey depth, in accordance with the selec- 
tion function (eqn. 19). Structure is most prominent in the 
first shell, at d sur = [30 — 55] ft, _1 Mpc (central left sphere). 
The structure contained in this shell also shows a clear af- 
filiation with the matter distribution in the local NBG 
volume. The compact massive concentration at I 220° 
is clearly connected to a dense region in the local "plane" . 
A superficial inspection of the angular galaxy distribution 
reveals the presence of strong dipolar and quadrupolar 
components, effecting considerable tidal forces. Note that 
both external shells display a rather strong concentration 
of galaxies in their southern hemisphere, in the vicinity of 
I w 180 — 200°. Similar but weaker contributions can also 
be recognized from the galaxy distribution in the shell be- 



tween d sur « 70 — 85/i _1 Mpc. Beyond d sur > 85/i _1 Mpc, 
however, the angular pattern appear to be considerably 
less pronounced. This is in line with the earlier findings 
that there were hardly noticeable tidal contributions from 
large distances. 

To see to what extent the depicted galaxy distribution 
can indeed be held responsible for most of the inferred 
tidal bulk flow and tidal shear, we have determined the 
corresponding bulk force gbuik (eqn. 4) and tidal shear %j 
(eqn. 5) evoked by the external galaxy distribution (r > 
30/i _1 Mpc). Since we do not have a continuous density 
field but the positions of a finite number of objects in our 
galaxy flux-limited and full mass distribution catalogs, the 
bulk acceleration on the LG is computed from the discrete 
equivalent. For a sample of galaxies at locations Xj, with 
an average number density n of selected objects, this leads 
to 



Qbulk 



Hf(n,A) 



Ann 



E 



1p{Xk) 



Xfc 

|xk| J 



(36) 



where ip(xk) is the sample selection function at dis- 
tance Xk, whose inverse functions as weighting factor. 
For practical reasons, comparison with the inferred bulk 
flow Vbuifc, we have translated the bulk acceleration into 
equivalent velocity units by means of the transformation 
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Fig. 18. Aitoff projections of the galaxy distribution of a NBG + PSCz rCDM catalog. The top-most panel shows the 
corresponding NBG distribution, d sur = [0, 30]/i -1 Mpc. The subsequent 4 panels depict the external galaxy distribution enclosed 
by the shells defined by survey depth d SU r = [30 — 55]/i _ Mpc, d sur = [55 — 70]ft _1 Mpc dsur- = [70 — 85]/i _1 Mpc and dsur = 
[85 — 100] h~ Mpc. Bottom panel: convergence gravity dipole and quadrupole. By means of symbols we have indicated the track 
of the gravity dipole and quadrupole eigenvector directions (cf. eqn. 36 and eqn. 37) on the sky, by radially expanding outward 
the survey depth dsur in steps of l/i -1 Mpc, from dsur = 30/i _1 Mpc to d SU r = 100/i _1 Mpc. For comparison the same symbol, but 
then enlarged, indicates the directions for the corresponding bulk velocity and velocity shear eigenvector directions. Diamond: 
dipole. Triangle: stretching component shear. Square: middle component shear. Star: compressional component shear. 
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Fig. 19. Cumulative gravity dipole and tidal gravity shear. As a function of survey depth d aur the Cartesian components of the 
gravity dipole, g x , g y and g z is followed (in equivalent velocity unit km s _1 ). The symbols at d sur indicate the corresponding 
velocity dipole. Lefthand panel: ACDM. Righthand panel: rCDM. 



Hf(n,A)/^ilH 2 . The equivalent "discrete" expression for 
the external tidal shear is 



Hf(n,A) 



E 



1 3 Xfci %kj 



(37) 



For the ACDM mock galaxy sample depicted in Fig 16 
we determine the gravity dipole by computing for a set 
of spherical external shells the resulting bulk flow ac- 
celeration (eq. 36) and the gravity quadrupole by com- 
puting external tidal shear (eq. 37). Recently, a simi- 
lar approach was followed by Teodoro et al. 2004. The 
spherical shell volumes are denned by an inner radius 
Tinn — 30/i _1 Mpc and an outer radius r out . The width 
of the shell is gradually enlarged by increasing r out from 
r ut = 30/i _1 Mpc to r out — 100ft. _1 Mpc. The convergence 
of the resulting gravity dipole direction on the sky can 
be observed in Fig. 17. The small red diamonds are con- 
sistently located near I m 230°, and converge at a sky 
location close to the direction of the velocity dipole (large 
red diamond). To get an idea of the amplitudes involved, 
Fig. 19 (top panel) shows the development of the cumula- 
tive gravity dipole as a function of external distance d sur 
(30/i~ 1 Mpc < dsur < 100/i _1 Mpc). By means of symbols 
the corresponding velocity dipole, for each of the three 
directions x, y and z, are inserted at the outer radius of 
d S ur ~ 100/i -1 Mpc. Note that we have restricted ourselves 
to a case study example. A more extensive and proper as- 
sessment, including a proper error estimate of both gravity 
dipole and quadrupole as well as the bulk and shear flow, 
is beyond the scope of the present argument. This issue, 
involving the shot noise influence on gravity dipole and 
quadrupole and the role of FAM uncertainties on the ve- 
locity flow components, will be treated in more detail in 
a forthcoming study. 

We see that in the x-, y and z-directions of the gravity 
and velocity dipoles are in reasonable agreement, within a 



margin of w 30 km s _1 . This observation justifies our ex- 
pectation that the dipole can be estimated to sufficient ac- 
curacy from the surrounding external galaxy distribution. 
The dipole may then be estimated from the surrounding 
external galaxy distribution, so that the latter can be in- 
voked to correct for the influence of the external tidal field 
in the dynamics of the local volume. 

The situation is comparable for the cumulative tidal 
shear, in terms of its three eigenvalues and eigenvectors. 
Also the gravity quadrupole appears to converge relatively 
smoothly towards the velocity shear. This may be inferred 
from the plotted directions of the eigenvectors en, £7-2 
and eq-3 of the tidal shear %j. They are indicated by 
means of three symbols, the triangle corresponding to the 
stretching component T\, the star the middle component 
T2 and the square the compressional component 73. The 
tidal shear tensor wanders extensively across the "sky" as 
we push the outer radius of the external shell outward, 
as is shown by the paths of the corresponding eigenvec- 
tors. Interestingly, once the shell radius starts to approach 
100/i _1 Mpc, each of the eigenvectors appear to converge 
near the location of the corresponding stretching, cen- 
tral and compressing velocity shear tensor eigenvectors. 
However, also here we notice significant deviations in in- 
dividual cases. 

For comparison, we can appreciate the role of the ex- 
ternal tidal field on local dynamics for the case of the 
rCDM cosmology. Figure 18 combines the galaxy sky dis- 
tribution for a tCDM mock galaxy sample, in the same 
radial shells as in Fig. 16. The final frame shows the Aitoff 
projection of the gravity dipole and gravity quadrupole 
eigenvectors for a set of gradually increasing radial shells. 
From the galaxy sky distribution in the four slices at 
sampling depths d sur = [0 — 30], [30 — 55], [55 — 70] and 
[85 — 100] ft, Mpc we notice that these involve consider- 
ably more isotropic distributions. Hardly any prominent 
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patterns can be discerned in the sky distribution. This is 
expressed in a more erratic wandering of gravity dipole 
and quadrupole directions (lower frame, fig. 16). This also 
implies a more substantial contribution of shot noise ef- 
fects. The latter represent a major source for deviations 
between the velocity dipole and shear flow quadrupole 
and the gravity dipole and quadrupole. The smaller co- 
herence length of the tCDM fluctuations and the more 
randomly oriented contributions by the individual exter- 
nal matter concentrations may therefore be directly re- 
lated to the lower level of coincidence between velocity and 
gravity directions than in the case of the more prominent 
anisotropies in the ACDM cosmology. In this, we have to 
realize that the amplitude of dipole and quadrupolar con- 
tributions between the two scenarios are not too different 
(cf. table 5). The less prominent anisotropies in the tCDM 
catalogues are therefore compensated by a higher average 
matter density. 

8.8. Multipole Components: Summary 

The above results reassure the fact that the external tidal 
field can be well characterized by its main multipole com- 
ponents, the bulk flow and velocity shear. In terms of mul- 
tipole amplitude convergence, these results show a better 
agreement for the ACDM model than for the tCDM one. 
This is due to the intrinsic characteristics of both cos- 
mic models. As has been discussed in Section 4, and may 
be directly appreciated from Fig. 2, the relatively lower 
amplitude of the rCDM perturbations is compensated 
by a higher mass content. It leads to an equally strong 
external gravitational influence. On the other hand, the 
smaller spatial coherence of density features in the rCDM 
scenario causes the orientation of the gravity dipole and 
quadrupoles to be rather jittery. The direction of the cu- 
mulative gravitational force in the rCDM scenario wan- 
ders erratically over the sky as we move further out from 
the local volume. This differs from the situation in the 
ACDM samples, where we observe a consistent, system- 
atic and coherent convergence towards the final dipole di- 
rection. 

The above results confirm the fact that the external 
tidal velocity field can be well characterized by its main 
multipole components, the bulk flow and velocity shear. 
This depends to some extent on the cosmology. In terms of 
multipole amplitude convergence, these results show a bet- 
ter agreement for the ACDM than for the rCDM model. 

9. Conclusions 

In this work we have applied the FAM technique to con- 
struct model velocity fields using mock catalogs resem- 
bling the NBG and IRAS-PSCz galaxy catalogs. The mock 
catalogs were extracted from ./V-body simulations in which 
the central observer mimics some of the properties of the 
Local Group environment. Comparing FAM velocities ob- 
tained from the NBG mock catalogs with those obtained 
from the larger PSCz mock catalogs and, finally, to the N- 



body velocities, allowed us to quantify the importance of 
the gravity field generated by the mass distribution within 
and beyond the LS. 

Neglecting the mass distribution outside the LS leads 
to a systematic underestimate of the gravity field. The 
amplitude of this bias depends on the amount of power 
on scales larger than the LS, and thus on the cosmological 
models. In a tCDM universe model peculiar velocities are 
<~ 20% smaller than the true ones. In the case of a ACDM 
model, which has more power on large scales, model ve- 
locities are underestimated by <~ 35%. 

The results of the described FAM analyses are encour- 
aging in the sense that the presently available all-sky, flux 
limited catalogs such as PSCz appear to be capable of ac- 
counting for the major share of the velocity field on the 
scale of the Local Supercluster. While the 30/i _1 Mpc re- 
stricted NBG sample showed a substantial deficiency in 
its capacity to generate the local cosmic motions, in par- 
ticular in the case of the ACDM Universe models, in both 
cases the 100/i~ 1 Mpc mock samples appear to embody 
nearly all matter concentrations responsible for the gener- 
ated velocities in our local (NBG catalog) neighbourhood. 

Also we notice a telling difference between the per- 
formance of both FAM30 and FAM100 reconstructions for 
the case of the ACDM cosmology catalogs on the one hand 
and the rCDM model catalogs on the other hand. The fact 
that the ACDM model involves substantially more power 
on large scales, r > 30/i~ 1 Mpc, than the rCDM model is 
reflected in the better quality of the FAM 30 reconstruc- 
tions for the rCDM catalogs. The presence of substantial 
mass inhomogeneities with a scale in excess of that of the 
local Universe regions implies a larger external contribu- 
tion to the local velocity field. This is also borne out by 
the fact that for the ACDM catalogs we see a consider- 
able improvement in velocity field reconstruction quality 
going from the FAM30 to the FAM100 reconstructions (see 
Table 4), while this is far less so for the rCDM catalogs. 

Of course, whether the resulting models do indeed form 
an unbiased representation of the actual velocity field will 
to some extent also depend on whether the galaxy distri- 
bution in the flux limited galaxy catalogs does represent 
an unbiased reflection of the actual (external) mass dis- 
tribution surrounding the Local Supercluster resembling 
region. The results of recent studies (Verde et al. 2002, 
Lahav et al. 2002, Tegmark, Zaldarriaga & Hamilton 2001 
and Branchini, Dckcl & Sigad 2002) are quite encouraging 
in this respect. They seem to indicate, certainly on scales 
larger than « 5/i _1 Mpc, that both IRAS and 2dF galax- 
ies trace the underlying mass distribution in an unbiased 
fashion. 

Nonetheless, observations along the lines of the pre- 
sented mock catalogs seem to suggest that a proper anal- 
ysis of Local Universe dynamics based on a combination 
of information of local small-scale (peculiar) galaxy veloc- 
ities and a rough yet well-founded idea of the matter dis- 
tribution on scales of a few hundred 100ft-~ 1 Mpc may help 
us towards acquiring far more insight into the dynamical 
history of the emergence and assembly of the striking non- 



44 



E. Romano-Diaz et al.: External Tidal Impact on Local Supercluster Dynamics 



linear patterns we have discovered in the large scale matter 
distribution. Moreover, we have uncovered evidence that 
a meticulous point-to-point analysis of such velocity sam- 
ples may help towards modelling the total local force field, 
including a proper model for the external forces. 

When modeling the peculiar velocity of a LS look- 
alike region by only considering the matter distribution 
within 30/i _1 Mpc, the end product is a biased velocity 
field lacking of any large scale signature. This bias can 
be eliminated by accounting for the mass distribution be- 
yond the LS. Our experiments demonstrate that sampling 
the mass distribution out to scales of 100/i _1 Mpc, in a 
flux limited fashion, is sufficient to account for the large 
scale contribution to the peculiar velocities in our cosmo- 
logical neighborhood. More precisely, we have found that 
the cosmic velocity field within the LS, modeled by FAM 
using the mass distribution traced by PSOz galaxies out 
to 100ft. _1 Mpc is unbiased. The differences between true 
and FAM velocity field are random and mainly occur in 
high density environments which are dominated by virial 
motions that are not modeled correctly by FAM. 

The gravity and velocity fields generated by the mass 
distribution beyond the scale of the Local Supercluster 
are well characterized by their bulk flow and shear com- 
ponents. Therefore, one can obtain an unbiased model ve- 
locity field by superimposing a local model velocity field 
within the Local Supercluster to the bulk flow and shear 
components of the velocity field generated by the mass 
distribution between 30 and 100/i _1 Mpc. 

These considerations suggest that velocity models 
which only consider the dynamics within the Local 
Superclusters might have been affected by systematic er- 
rors. In particular, our work suggests that, when com- 
pared with observed velocities, they might have underesti- 
mated the value of the density parameter, Q m , by 15-25 %. 
However, the analysis of Shaya, Tully and Peebles (1995), 
based on the galaxy distribution in the Local Supercluster, 
shows that a lower, not a larger, value of fl m is found when 
complementing the local mass distribution with the large 
scale one traced by rich Abell clusters. A more precise 
evaluation of this bias will be performed in a future work 
in which we will perform the same analysis presented here 
using a new set of mock catalogs that are constrain to re- 
produce the distribution of the mass in our local Universe 
(see e.g. van de Weygaert & Hoffman 1999, Mathis et al. 
2002, Klypin et al. 2003). 

Furthermore, our analysis shows that all model veloc- 
ity fields of the Local Supercluster which are based on 
the PSCz catalog (e.g. Branchini et al. 1999, Schmoldt 
et al. 1999, Valentine et al. 2000, Sharpe et al. 2001) are 
free from systematic biases arising from having neglected 
the large scale contribution from scales beyond its realm. 
Moreover, since the IRAS PSCz survey is considerably 
deeper than 100/i _1 Mpc, it is reasonable to assume that 
the PSCz catalog can be used to predict unbiased veloci- 
ties well beyond our Local Supercluster that, if compared 
with observed galaxy peculiar velocities, can discriminate 



among different cosmologies characterized by different val- 
ues of Q m , like the ACDM and rCDM models. 

The plausibility of this hypothesis has been recently 
confirmed by the analysis of Hoffman et al. (2001) that 
shows that the bulk and shear components of the external 
velocity field in the local universe inferred from the pecu- 
liar velocities in the Mark III catalog (Willick et al. 1997a, 
1997b) are qualitatively consistent with those expected 
from the mass distribution traced by IRAS PSCz galaxies. 
On the other hand, the claim on the basis of the SMAC 
cluster peculiar velocity sample (Hudson et al. 2004) of 
an extra 225 km s _1 bulk flow component generated by 
matter concentrations on a scale exceeding 100/i _1 Mpc 
should issue some caution with respect to claims of having 
accounted for all external influences on the local cosmic 
flow. 

Coupling the local velocity model provided by FAM 
to the large scale contribution provided by linear theory 
allows to obtain a model velocity field which is unbiased, 
nonlinear and fast to compute. This means that, for the 
first time, we are in the position of performing a large 
number of experiments aimed at studying the nonlinear 
evolution of cosmic structures, such as filaments and clus- 
ters, and explore the role of tidal fields during their gravi- 
tational collapse. This relates to the observation that fila- 
ments are forming as a consequence of anisotropic collapse, 
induced a comprcssional tidal force acting perpendicular 
to the "axis" of the filament. By tracing out the coherent 
paths of the compressional modes of the primordial tidal 
field one can identify the sites of the later nonlinear fila- 
ments (Bond et al. 1996, see Van de Weygaert 2002). In 
turn this is directly related to cluster locations: the strong 
primordial tidal shear is the result of a local quadrupolar 
mass distribution. The corresponding overdensities tend 
to evolve into rich clusters, explaining the intimate link of 
clusters and filaments in the cosmic web. 

Finally, it is worth stressing that in this work we 
have neglected the fact that we measure galaxy redshifts 
rather then positions. By means of an elegant formal- 
ism, Phelps 2000 demonstrated the feasibility of work- 
ing out the action principle in redshift space. With re- 
spect to FAM, Nusser & Branchini (2000) have shown 
that it can be easily implemented in redshift space and 
Branchini, Eldar & Nusser (2002) demonstrated that it 
performs equally well in real and redshift space. Therefore, 
our unbiased, nonlinear model velocity field also allows to 
perform an accurate correction for redshift space distor- 
tions and thus lead to a precise reconstruction of the mass 
distribution in real space. Mapping the mass in the local 
universe down to nonlinear scales and comparing it with 
the distribution of baryonic mass (in form of stars or dif- 
fuse, ionized gas) is of considerable astrophysical interests 
as it will constrain and help understanding the process of 
galaxy formation and evolution within the Universe 
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